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ABSTRACT 


A  new  technique  for  imaging  the  2D  transport  of  free  charge  in 
semiconductor  structures  is  used  to  directly  map  electric  field  distributions  in 
operating  devices.  Direct  transport  imaging  is  demonstrated  in  a  scanning 
electron  microscope,  using  an  optical  microscope  and  a  high  sensitivity  charge 
coupled  device.  Transport  behavior  under  the  combined  influence  of  both 
diffusion  and  drift  is  predicted  by  modeling  the  drift  and  diffusion  in  2D  following 
generation  at  a  point  source.  This  is  the  first  demonstration  of  a  technique  that 
allows  the  mapping  of  the  electric  field  by  determining  not  only  the  direction  but 
especially  the  magnitude  of  the  electric  field  with  high  resolution.  The  measured 
results  show  excellent  agreement  with  theoretical  predictions  simulated  with 
COMSOL  software. 

The  transport  imaging  technique  also  allows  measurement  of  the  contact 
resistance  in  a  new  way  that  is  non-destructive  and  based  on  a  two-point  contact 
only.  The  technique  illustrates  the  device’s  characteristics  by  determining  the 
exact  activation  point  of  the  diode  and  the  deviations  from  an  ideal  l-V  behavior. 
The  method  is  extremely  useful  since  the  complexity  and  miniaturization  of 
current  devices  do  not  allow  for  multiple  wiring  that  standard  four  point 
measurement  demands. 

Finally,  a  suggestion  for  further  research  of  the  effects  of  electromigration 
by  using  the  direct  transport  imaging  technique  is  offered.  The  latter  is  a  subject 
of  high  importance  in  electronic  device  reliability. 
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I.  INTRODUCTION 


A.  DIRECT  TRANSPORT  IMAGING  AND  ELECTRIC  FIELD 

VISUALIZATION  IN  SEMICONDUCTOR  PLANAR  STRUCTURES 

Electro-optics  is  the  field  of  systems  that  convert  photons  to  electrons 
regardless  of  the  wavelength  [1],  Detectors  are  used  across  the  full  range  of 
electromagnetic  (EM)  spectrum,  with  special  military  interest  in  the  near  or  short 
wave  infrared  [(SWIR),  0.7-2  pm],  in  the  midwave  infrared  [(MWIR),  3-5  pm]  and 
in  the  long-wavelength  infrared  [(LWIR),  8-12  pm]  regions.  The  sensitivity  of 
these  detectors  of  photons  and  of  high  energy  particles  depend  on  the 
generation  and  subsequent  motion  and  collection  of  free  charge.  Similarly,  a 
wide  range  of  devices  depend  upon  the  fabrication  of  contacts  that  can  uniformly 
inject  and  remove  charge.  New  efforts  in  materials  and  devices  must  be 
accompanied  by  improved  understanding  of  local  transport  behavior  [2]. 

In  particular  it  is  of  high  importance  to  be  able  to  directly  image  the 
injection  and  motion  of  free  charge  with  high  resolution  since  a  macroscopic 
current-voltage  measurement  is  based  on  the  averaging  of  local  transport 
behavior.  Various  techniques  have  been  adopted  in  order  to  study  the 
luminescent  phenomena  and  transport  properties  of  structures  used  in  sensor 
device  manufacture.  The  work  described  in  this  thesis  is  based  on 
cathodoluminescence  scanning  electron  microscopy,  which  allows  for  spatial 
characterization  of  light  emission. 

In  this  thesis,  a  technique  for  imaging  the  2D  transport  of  free  charge  in 
planar  semiconductor  structures  has  been  developed  in  our  laboratory  [2].  This 
technique  is  used  to  image  the  motion  of  minority  carriers  in  a  heavily  doped 
p-type  AIGaAs/GaAs/AIGaAs  heterostructure,  under  a  range  of  applied  electric 
fields.  The  laboratory’s  scanning  electron  microscope  (SEM)  is  a  JEOL  840A  with 
an  internal  optical  microscope,  connected  to  an  external  thermoelectrically 
cooled  charge  coupled  device  (CCD)  camera.  An  electron  beam  is  generated 
through  thermionic  emission  and  interacts  with  the  specimen,  creating  electron- 
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hole  pairs  and  finally  the  emission  of  photons.  By  generating  charge  at  a  fixed 
point,  it  is  possible  to  monitor  the  subsequent  drift  and  diffusion  by  imaging  the 
distribution  of  the  free  carrier  recombination  path  [2], 

The  main  objective  of  this  thesis  is  to  demonstrate  electric  field  mapping 
by  determining  not  only  the  direction  but  especially  the  magnitude  of  the  field 
locally.  In  order  to  determine  the  extent  to  which  transport  imaging  can  be  used 
to  map  local  electric  fields  with  high  resolution,  charge  motion  as  a  function  of 
applied  electric  field  has  been  imaged.  Transport  behavior  under  the  combined 
influence  of  both  diffusion  and  drift  is  predicted  by  modeling  the  drift  and  diffusion 
in  2D  following  generation  at  a  point  source.  Subsequent  measurements  of 
electric  field  profiles  have  been  done  in  both  uniform  (parallel  plate  structure)  and 
non-uniform  regions. 

With  the  method  mentioned  above  it  is  also  possible  to  measure  the 
contact  resistance  in  a  new  way.  In  order  to  be  able  to  define  the  applied  electric 
field  strength  in  an  area  where  drift  measurements  take  place,  a  four  point 
measurement  is  an  efficient  way.  The  disadvantage  of  the  method  is  that  the 
complexity  and  small  size  of  newly  manufactured  semiconductor  planar  devices 
does  not  allow  multiple  wiring  for  additional  current  carrying  and  voltage 
measuring  probes.  By  applying  the  electron  beam  in  the  uniform  region  between 
two  metal  contacts  superjacent  to  semiconductor  material,  study  of  the  transport 
behavior  under  the  influence  of  both  diffusion  and  drift  at  the  point  source  in  a 
range  of  applied  voltage  allows  accurate  measurement  of  the  contact  resistance 
with  only  two-point  wiring.  The  technique  illustrates  the  device’s  characteristics 
by  determining  the  “turn  on  voltage”  of  the  diode  and  further  deviations  from  an 
ideal  linear  l-V  behavior. 

Finally,  an  investigation  of  electric  field  distributions  relevant  to 
electromigration  has  been  started  for  potential  future  applications.  In  particular, 
an  effort  to  generate  electromigration  in  a  specimen  has  been  performed  in  order 
to  study  how  the  electric  field  profiles  of  specific  devices  are  affected  in  an 
Electronic  Warfare  environment  (e.g.,  radiation  from  electronic  attack). 
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Nevertheless,  the  characteristics  of  the  interaction  between  current  carriers  and 
migrating  atoms  are  not  entirely  understood.  Further  applications  of  the  direct 
transport  imaging  technique  can  be  used  for  future  study  of  the  electromigration. 
The  latter  is  currently  a  subject  of  high  importance  in  electronic  device  reliability. 

B.  MILITARY  RELEVANCE 

The  detector  is  the  component  of  the  sensor  system  which  transforms  the 
optical  signal  into  an  electrical  signal.  The  spectral  response  is  determined  by  the 
detector  characteristics  and  the  operating  temperature.  The  sensitivity  of  the 
material  used  (threshold  wavelength)  is  determined  by  its  energy  bandgap  (Eg). 
The  wavelength  of  a  photon  has  to  be  short  enough  so  that  the  photon  energy 
exceeds  a  material’s  bandgap  energy  level  for  absorption.  Intrinsic  detectors  are 
constructed  with  semiconductor  crystals  with  no  introduced  impurities  or  lattice 
defects,  while  extrinsic  detectors  have  impurities  introduced  in  order  to  change 
the  effective  bandgap  energy  [1]. 

In  the  case  of  photon  or  quantum  detectors,  the  change  in  charge  carrier 
state  changes  the  electrical  properties  of  the  material  indicating  the  amount  of 
incident  optical  power.  Depending  on  their  characteristics,  photon  detectors  can 
be  distinguished  as  photoconductive,  photovoltaic  and  photoemissive  devices  [1]. 
From  all  the  above  it  is  obvious  that  the  key  component  of  each  type  of  detector 
is  the  semiconductor  material  that  converts  the  absorbed  photon  in  the  desired 
electrical  signal. 

Military  systems  for  target  detection,  ranging  and  remote  observation  are 
highly  dependent  upon  advanced  photonic  and  optoelectronic  devices.  It  is 
critical  that  these  devices  are  highly  sensitive,  reliable  and  cost  effective  in 
operation.  The  military,  along  with  all  other  sectors  of  the  economy,  has  benefited 
greatly  from  the  advances  in  semiconductor  technology  over  recent  decades. 

New  demands  constantly  require  new  devices.  Recent  examples  would 
include  the  production  of  lasers  and  LEDs  for  the  blue  and  UV  parts  of  the 
spectrum,  improved  radiation  sensors  for  room  temperature  operation,  high 
temperature  materials,  such  as  SiC  and  GaN,  for  demanding  environments  and 
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the  renewed  interest  in  terahertz  detectors  for  explosive  detection.  In  all  these 
cases,  new  materials  and  new  contact  technologies  must  be  developed. 

In  this  development  process,  characterization  tools  are  required  that 
provide  the  needed  information  in  a  cost  effective,  timely  manner.  The 
characterization  tool  that  is  developed  and  tested  in  this  thesis  will  provide  a 
means  for  direct  characterization  of  new  nanostructures  and  other  devices  that 
are  being  considered  for  next  generation  electronics.  In  addition,  because  it 
utilizes  a  scanning  electron  microscope,  which  is  an  existing  tool  for  high 
throughput,  commercial  applications,  the  technique  has  the  potential  to  provide 
rapid,  large  scale  transport  characterization.  This  could  reduce  the  cost  of  new 
devices  by  providing  early,  device-critical  feedback  with  a  characterization  tool 
that  requires  minimal  device  processing. 

Finally,  device  reliability  in  operation  requires  a  full  understanding  of 
material  defects  as  well  as  the  interaction  of  material  with  the  environment,  such 
as  in  the  radiation  behavior  of  devices  in  space.  This  technique  has  the  potential 
to  provide  detailed  defect  characterization  of  devices  during  operation.  By 
mapping  the  electric  field  on  a  micron  scale,  early  indication  can  be  obtained  of 
future  failure  modes.  As  will  be  demonstrated,  this  technique  has  the  potential  to 
measure  local  electric  fields  at  about  10%  variation,  providing  a  very  sensitive 
measure  of  developing  defect  behavior  that  could  ultimately  affect  device  scale 
performance.  This  capability  is  essential  due  to  miniaturization  of  the  current 
devices  and  the  complexity  of  the  influence  that  various  defects  have  on  their 
electrical  behavior. 

C.  THESIS  OVERVIEW 

Chapter  I  states  the  purpose  of  the  thesis  providing  an  overview  of  the 
experimental  technique  of  direct  transport  imaging  in  order  to  locally  map  and 
measure  the  electric  field  in  semiconductor  planar  structures.  Chapter  II  provides 
background  information  concerning  the  physics  of  semiconductor  materials  and 
luminescent  phenomena  in  semiconductors  and  an  overview  of 
cathodoluminescence.  Chapter  III  describes  the  experimental  approach  of  the 
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cathodoluminescence  scanning  electron  microscopy  followed  in  the  current  work. 
It  includes  a  description  of  the  laboratory’s  equipment,  the  specimen  used,  the 
overall  procedure  and  the  mathematical  models  used.  In  Chapter  IV,  a  number  of 
cases  are  illustrated  in  order  to  demonstrate  the  validity  and  potential  of  the 
method  described  in  Chapter  III.  In  Chapter  V,  a  potential  method  of  evaluating 
the  contact  resistance  with  a  new  two-point  measurement  has  been  presented. 
Additionally,  in  the  same  chapter,  a  method  of  evaluating  the  electric  field 
magnitude  in  semiconductor  planar  structures  with  constricted  geometry  has 
been  presented.  Chapter  VI  refers  to  electromigration  issues  and  provides 
conclusions  and  suggestions  for  future  work. 
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II.  BACKGROUND 


A.  CONDUCTION  IN  SEMICONDUCTORS 

In  semiconductors,  there  are  two  types  of  free  charge  carriers:  electrons 
and  holes.  A  hole  is  effectively  an  empty  electronic  state  in  the  valence  band 
(VB)  that  behaves  as  if  it  were  a  positively  charged  particle  free  to  respond  to  an 
applied  electric  field.  Free  electrons  that  exist  in  the  conduction  band  (CB)  can 
also  move  in  the  crystal  and  contribute  to  the  electrical  conduction  when  an 
electric  field  is  applied. 

When  small  amounts  of  impurities  are  introduced  into  an  otherwise  pure 
crystal  structure,  it  is  possible  to  obtain  a  semiconductor  in  which  the 
concentration  of  carriers  of  one  polarity  is  much  in  excess  of  the  other  type.  This 
kind  of  semiconductor  is  called  an  extrinsic  semiconductor.  If  the  impurities  have 
a  smaller  valency  than  the  equivalent  of  the  pure  semiconductor  material,  the 
result  is  an  excess  of  holes  over  electrons.  This  is  a  p-type  semiconductor.  A 
characteristic  example  is  the  p-type  material  created  by  Si  doped  with  B,  Si:B. 

Si  has  Z=14  and  electronic  structure  [Ne]3s^p^  and  B  has  Z=5  and 
electronic  structure  [He]2s^p\  It  can  be  seen  that  B  has  three  valence  electrons 
and  Si  four  valence  electrons.  So,  when  B  shares  its  valence  electrons  with  four 
Si  atoms,  one  of  the  bonds  has  a  missing  electron  which  is  a  hole.  This  effective 
empty  state  in  the  valence  band  results  in  the  characterization  of  B  as  an 
acceptor  (electron  acceptor)  impurity  as  it  accepts  an  additional  electron  in  order 
to  form  four  covalent  bonds  around  it.  The  primary  material  used  in  this  thesis  is 
GaAs:Be  which  is,  according  to  the  above,  a  p-type  semiconductor  material. 

When  an  electron  in  the  CB  meets  a  hole  in  the  VB,  the  electron  occupies 
the  hole  since  it  found  a  state  of  lower  energy.  This  recombination  process 
results  in  the  annihilation  of  an  electron  in  the  CB  and  a  hole  in  the  VB.  The 
excess  energy  of  the  electron  returning  from  CB  to  VB  in  certain  semiconductors 
such  as  GaAs  is  emitted  as  a  photon  [3].  It  is  this  emission  process  that  is 
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monitored  in  this  thesis  in  order  to  image  the  transport  of  free  charge  in  planar 
semiconductor  structures  under  the  influence  of  applied  electric  field. 

B.  DOUBLE-HETEROSTRUCTURE  (DH)  DEVICE 

A  junction  between  two  different  bandgap  semiconductors  is  called  a 
heterojunction.  A  semiconductor  device  structure  that  has  junctions  between 
different  bandgap  materials  is  called  a  heterostructure  device  [3]. 

GaAs  and  AIGaAs  are  two  semiconductors  with  the  characteristics  in 
Table  1. 


Semiconductor  Active  Layer 

Structure 

Bandgap 

GaAs 

DH 

Direct  (Eg»1 .4eV) 

AlxGai-xAs  (0<x<0.4) 

DH 

Direct  (Eg«2eV) 

DH=Double  Heterostructure 

Table  1.  Selected  Semiconductor  Materials. 


Element 

Symbol 

Z 

Electronic  Structure 

Gallium 

Ga 

31 

[Ar]3d^°4s^p^ 

Arsenic 

As 

33 

[Ar]3d^°4s^p^ 

Aluminum 

Al 

13 

[Ne]3s^p^ 

Boron 

B 

5 

[He]2s^p^ 

Beryllium 

Be 

4 

[He]2s^ 

Table  2.  Selected  Elements  of  the  Periodic  Table. 


Figure  1.  Double-Heterostructure  Diode.  {After  [3]}. 
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The  double  heterostructure  diode  in  Figure  1  has  an  n'^p  heterojunction 
between  n""  -AIGaAs  and  p-  GaAs  and  another  heterojunction  between  p-  GaAs 
and  p-  AIGaAs. 


Figure  2.  Energy  Band  Diagram  with  no  bias.  {After  [3]}. 


If  no  bias  is  applied  in  the  above  double  heterostructure,  the  energy  band 
diagram  can  be  seen  in  Figure  2.  The  Fermi  energy  (Ef)  level  is  uniform.  There  is 
a  potential  energy  barrier  eVo  for  electrons  in  the  CB  of  n"^  -AIGaAs  against 
diffusion  into  p-  GaAs.  AEc  is  a  potential  energy  barrier  that  prevents  any 
electrons  in  the  CB  in  p-  GaAs  moving  to  the  CB  of  p-  AIGaAs. 


um 


fonvard 


Figure  3.  Energy  Band  Diagram  with  forward  bias.  {After  [3]}. 
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When  a  forward  bias  is  applied,  the  potential  energy  barrier  eVois  reduced 
and  electrons  in  the  CB  of  r\^  -AIGaAs  are  injected  into  p-  GaAs  as  shown  in 
Figure  3.  Nevertheless,  the  energy  barrier  AEc  confine  these  electrons  in  the  p- 
GaAs  region  not  allowing  them  to  diffuse  to  the  p-  AIGaAs  region.  The 
recombination  of  injected  electrons  to  the  holes  that  already  exist  in  this  p-  GaAs 
layer  results  in  spontaneous  photon  emission. 

The  example  mentioned  above  is  a  characteristic  example  of  a  double 
heterostructure  (DH)  laser  diode  based  on  two  junctions  between  different 
semiconductor  materials  with  different  energy  bandgaps  (n-p-p  type  material). 
The  reason  for  describing  the  principles  of  operation  of  this  kind  of  diode  is  that 
the  experimental  approach  developed  in  this  thesis  is  based  on  a  solid  source 
molecular  beam  epitaxy-grown  AI0.90Ga0.10As/GaAs/AI0.90Ga0.10As  double 
heterostructure  sample,  so  it  is  important  to  be  able  to  understand  the  carrier 
concentration  in  the  GaAs  layer  in  order  to  study  the  transport  behavior  of  free 
charge. 

C.  CARRIER  TRANSPORT  PHENOMENA 

When  an  electric  field  Ex  is  applied  to  a  semiconductor  material,  electrons 
in  the  CB  move  opposite  to  the  direction  of  the  electric  field,  gain  energy  from  the 
field  and  lose  energy  from  collisions  due  to  thermal  vibrations.  This  process  is 
called  drift  of  the  electrons  in  an  applied  field.  Similarly,  the  holes  in  the  VB  also 
drift  in  an  applied  field  Ex  but  now  the  drift  is  in  the  direction  of  the  electric  field. 
Since  both  electrons  and  holes  contribute  to  electrical  conduction,  current  density 
J  is  given  by: 

J  =  enVde  +  epVdk  (1) 

where  n  is  the  electron  concentration  in  the  CB,p  is  the  hole  concentration  in 
the  VB,  e  is  the  electron  charge  and  Vde  ,  Vdh  are  the  drift  velocities  of  electrons 
and  holes  respectively  [3]. 

The  drift  velocities  of  electrons  Vde  and  holes  Vdh  in  response  to  an  applied 
electric  field  Ex  are  given  by: 
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Vde  =  JJeEl 


(2) 


Vdh  =  /uhEx  (3) 

where  jue  and  //a  are  the  electron  and  hole  drift  mobilities  [3]. 

Drift  mobility  is  the  constant  of  proportionality  between  the  drift  velocity  Vd 
and  the  applied  electric  field  Ex .  The  drift  mobility  jue  of  the  electrons  and  the 
drift  mobility  juh  of  the  holes  are  given  by: 


_  etc 
=  T 
Me 


(4) 


ezh 

^ 


nih 


(5) 


where  le  and  ih  are  the  mean  free  times  between  scattering  events  and  me* 
and  mh  are  the  effective  electronic  mass  and  effective  mass  of  a  hole 
respectively  [3].  The  conductivity  of  a  semiconductor  cr  is  then  given  by  [3]: 

cr  =  en/Ue  +  ep/Uh  (6) 

The  diffusion  coefficient  is  a  measure  of  the  ease  with  which  the  diffusing 
charge  carriers  move  in  the  medium.  Diffusion  is  the  random  process  by  which 
particles  move  from  high  concentration  regions  to  low  concentration  regions. 
Similarly,  drift  mobility  is  a  measure  of  the  ease  with  which  the  charge  carriers 
move  in  the  medium  due  to  an  applied  electric  field  [3].  The  above  two  quantities 
are  related  through  the  Einstein  relation: 


De  kT  ^  kT 

- = - ^  De  =  jUe - 

(7) 

jUe  e  e 

Dh  kT  _  kT 

—  =  —  ^  Dh  =  jUh  — 

(8) 

juh  e  e 

where  De  and  Dh  are  drift  coefficients,  //e  and  juh  are  drift  mobilities  of 
electrons  and  holes  respectively,  k  is  the  Boltzmann  constant  and  T  is  the 
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absolute  temperature.  The  —  quantity  facilitates  the  use  of  units  of  eV  for  the 

e 

determination  of  kT . 

It  can  be  seen  from  equations  (7)  and  (8)  that  the  diffusion  coefficient  is 
proportional  to  the  temperature  and  mobility.  This  is  reasonable  since  when 
temperature  is  increased,  diffusion  will  be  enhanced  too. 

The  diffusion  coefficient  is  also  affected  by  the  distance  /  that  a  free 
electron  (or  hole)  covers  before  scattering,  in  time  t  by  [3]: 

D  =  -  (9) 

T 

Based  on  equation  (9),  the  diffusion  length  of  holes  is  given  by  [3]: 

Lh  =  ^iDh)x{Th)  (10) 

In  a  p-type  semiconductor,  holes  are  the  majority  carriers  and  electrons 
are  the  minority  carriers.  Minority  carrier  concentration  in  this  case  is  termed  as 
npo  where  the  subscript  “po”  refers  to  p-type  semiconductor  at  thermal 
equilibrium.  Thermal  equilibrium  means  that  the  mass  action  law  is  obeyed  [3]: 

tlpoppo  ~  fli  (11) 

where  m  is  the  intrinsic  carrier  concentration. 

It  should  be  mentioned  that  in  equation  (10),  the  term  Th  stands  for  the 
minority  carrier  lifetime  which  is  the  average  time  that  a  hole  exists  in  the  VB 
from  its  generation  to  its  recombination  with  an  electron.  Recombination  time  Th 
depends  on  the  semiconductor  material,  impurities,  crystal  defects,  temperatures 
and  so  forth,  and  there  is  no  typical  value  to  quote  (varies  from  nanoseconds  to 
seconds).  In  heavily  doped  semiconductor  material  Ts  ^  Trecom.  which  means  that 
the  scattering  time  is  almost  identical  to  the  recombination  time. 
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D.  LUMINESCENCE  PHENOMENA  IN  SEMICONDUCTORS 

Luminescence  is  the  process  of  photon  emission  when  a  solid  is  supplied 
with  a  certain  form  of  energy  [4],  Different  names  are  used  depending  on  the 
source  of  the  excitation.  In  photoluminescence  the  excitation  is  due  to  absorption 
of  photons.  In  cathodoluminescence  the  excitation  is  being  done  by  absorption  of 
energetic  electrons  or  cathode  rays.  In  chemiluminescence  the  required  energy 
for  the  excitation  is  supplied  by  chemical  reaction.  In  electroluminescence  the 
excitation  is  due  to  the  application  of  an  electric  field.  In  semiconductors, 
luminescence  is  generally  due  to  the  radiative  recombination  of  electron-hole 
pairs. 

The  emission  of  photons  in  luminescence  is  due  to  the  electronic  transition 
between  the  initial  energy  state  Ei  and  the  final  energy  state  Ef.  The  energy  and 
the  wavelength  of  the  emitted  photon  can  be  found  from  the  equation: 

E  =  Ef-Ei  =  hv  =  —  (12) 

X 

where  v  and  X  are  the  frequency  and  wavelength  of  the  transmitted  photon,  h 
is  Planck’s  constant  and  c  is  the  speed  of  light. 


E  E 


Figure  4.  Direct  and  Indirect  Energy  Band  Diagram.  {After  [3]}. 

In  the  first  case  illustrated  in  Figure  4,  GaAs  is  a  direct  bandgap 
semiconductor.  During  the  recombination  of  electron-hole  pair,  the  electron  drops 
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from  the  bottom  of  the  CB  to  the  top  of  the  VB  without  any  change  in  its  k  value 
{k  is  the  wave  vector)  and  the  momentum  is  conserved  (hk)  [3].  The  result  of  the 
recombination  is  the  emission  of  a  photon. 

In  the  second  case  illustrated  in  Figure  4,  Si  is  an  indirect  bandgap 
semiconductor  since  the  minimum  of  the  CB  is  not  directly  above  the  maximum 
of  the  VB.  If  there  is  an  electron  at  the  bottom  of  the  CB,  it  cannot  recombine 
directly  with  a  hole  at  the  top  of  the  VB  because  momentum  is  not  conserved 
{kcb  ^  kvh )  [3].  In  order  for  an  electron  to  recombine  with  a  hole  in  such  a  material, 
an  intermediate  recombination  center  with  energy  Er  should  be  used,  as  shown  in 
Figure  5  below. 


E 


Figure  5.  Indirect  Energy  Band  Diagram  with  Recombination  Center. 

{After  [3]}. 

In  this  case,  part  of  the  energy  of  the  electron  is  lost  by  the  emission  of 
phonons  that  correspond  to  lattice  vibrations.  Nevertheless,  in  some  indirect 
bandgap  semiconductors,  the  recombination  process  via  a  recombination  center 
results  in  photon  emission  [3].  Such  a  process  has  a  much  lower  probability  of 
occurrence  than  direct  transitions  and  so  fundamental  emission  in  indirect-gap 
semiconductors  is  relatively  weak. 


14 


E.  CATHODOLUMINESCENCE 

Electron  beam  excitation  may  produce  orders  of  magnitude  greater  carrier 
generation  rates  than  typical  optical  excitation.  The  generation  factor,  which  is 
the  number  of  electron-hole  pairs  generated  per  incident  electron  beam,  is  given 
by: 

G^Eb{\-y) 

Ei 

where  Eb  is  the  electron  beam  energy,  Ei  is  the  ionization  energy  (the  energy 
required  for  the  formation  of  the  electron-hole  pair)  and  y  is  the  fractional 
electron  beam  energy  loss  due  to  the  backscattered  electrons  [4], 

The  ionization  energy  Ei  is  related  with  the  band  gap  of  the  material  by: 

Ei  =  2.Mg+M  (14) 

In  equation  (14),  0<M  <^  eV.  M  depends  on  the  material  and  is  independent  of 

the  electron  beam  energy  [4], 

The  local  generation  rate  of  carriers  is: 

g{r,z)  =  {g)—  (15) 

e 

where  (g)  is  the  normalized  distribution  of  the  ionization  energy  in  the 

generation  volume.  It  is  the  electron  beam  current  and  e  is  the  electronic 
charge  [4], 

The  local  generation  rate  of  carriers  has  been  determined  experimentally 
for  silicon  and  a  universal  depth-dose  function  g(z)  is  given  by: 

g(z)  =  0.60  +  1.21z-12.40zV5.69z^  (16) 

Equation  (16)  represents  the  number  of  electron-hole  pairs  generated  per 
electron  of  energy  E  per  unit  depth  per  unit  time.  Graphically  it  can  be 
represented  as  shown  in  Figure  6  [4].  GaAs  will  show  a  similar  behavior  for 
varying  values  of  incident  beam  energy. 
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Figure  6.  The  depth-dose  curves  for  Si.  {From  [5]}. 

F.  CATHODOLUMINESCENCE  INTENSITY 

The  CL  intensity  represents  the  number  of  photons  emitted  per  unit  time 
and  is  given  by  [4]: 

LcL{r)  =  \f^^d^r  (17) 

V 


An(r^ 

where  r.,  is  the  radiative  recombination  lifetime,  — ^  is  the  radiative 

Trr 

recombination  rate,  An(r)  is  the  excess  minority  carrier  density  and  f  is  a 
function  containing  correction  parameters  of  the  CL  detection  system  and  factors 
that  account  for  the  fact  that  not  all  photons  generated  in  the  material  are  emitted 
(due  to  optical  absorption  and  reflection  losses)  [4], 

For  a  spherically  symmetric  distribution  far  from  the  point  source: 

r 

An{r)  =  C—  (18) 

r 


where  C  is  a  constant  and  L  =  4^  is  the  minority  carrier  diffusion  length  [4],  It 
should  be  noted  that  equation  (18)  refers  to  diffusion  in  three  dimensions.  In  this 
thesis  where  the  sample  is  a  thin  layer  and  a  2D  distribution  from  the  point 
source  is  used,  equation  (18)  becomes: 


An(r)  =  CKo 


kLj 


(19) 
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where  Ko  is  the  zeroth  order  modified  Bessel  function  of  the  second  kind  (the  2D 
case  is  described  in  Chapter  III  of  the  current  thesis). 

The  luminescence  intensity  due  to  the  radiative  recombination  in  a  layer  of 
thickness  dz  at  a  depth  z  is: 

z 

A«(r)  GIhLe  ^dz 

LcL{z)dz  oc - oc -  (20) 

Trr  eVrrD 

More  specifically  the  intensity  is  given  by: 

LcL=f44.—  (21) 

eXrr 

where  /b  is  a  constant  that  accounts  for  parameters  of  the  CL  detection  system 
such  as  the  overall  collection  efficiency  of  the  light  collector,  fA  and  fR  are 
factors  that  account  for  absorption  and  internal  reflection  losses  [4], 

An  understanding  of  the  cathodoluminescent  intensity  is  crucial  because 
the  work  presented  in  this  thesis  is  based  on  the  direct  transport  imaging  of  free 
charge  carriers  by  simultaneously  evaluating  the  luminescent  intensity  due  to 
radiative  recombination  in  a  point  source.  The  values  of  the  intensity  that  can  be 
acquired  by  using  an  appropriate  2D  mathematical  model  for  DC  and  AC  electric 
fields  will  be  used  in  order  to  be  able  to  determine  the  values  of  the  electric  field 
locally  in  planar  semiconductor  material. 
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III.  EXPERIMENTAL  APPROACH 


A.  DIRECT  TRANSPORT  IMAGING 

It  is  very  important  to  be  able  to  directly  image  the  injection  and  the  motion 
of  free  charge  in  planar  semiconductor  structures,  with  high  resolution.  By  doing 
this,  much  more  can  be  learned  than  from  macroscopic  current-voltage 
measurements  that  are  based  on  the  average  of  local  transport  behavior.  The 
small  size  of  current  semiconductor  devices  used  in  industry,  the  structure 
complexity,  the  role  of  the  material  non-uniformity,  the  effects  of  various  defects 
and  the  necessity  of  a  non-destructive  approach  are  some  of  the  reasons  that 
make  the  direct  transport  imaging  technique  an  efficient  way  to  study  the  electric 
behavior  of  semiconductor  devices,  particularly  sensors  and  detectors  that 
depend  on  the  collection  of  photo-generated  charge. 

The  direct  transport  imaging  technique  described  in  this  thesis  is  based  on 
the  capability  to  image  the  motion  of  free  carriers  in  planar  semiconductor 
structures  under  a  range  of  applied  electric  fields  via  the  luminescence  emission. 
This  is  illustrated  in  Figure  7  below. 


Figure  7.  Spot  mode  imaging  in  a  p-type  AIGaAs/GaAs/AIGaAs 

heterostructure  in  two  different  cases  of  applied  electric  field 
[(a)  and  (b)  216  pm  X  190  pm]. 
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As  in  CL,  the  recombination  of  electron  -  hole  pairs  created  by  the 
interaction  of  the  electron  beam  with  the  planar  structure  results  in  photon 
emission.  The  difference  with  the  conventional  CL  technique  is  that  the  electron 
beam  is  held  fixed  over  a  predefined  location  [2],  By  doing  this,  the  generation 
point  (of  e-h  pairs)  is  held  fixed  too.  The  largest  fraction  of  the  photon  signal  due 
to  the  recombination  is  created  at  or  near  the  point  of  charge  generation. 
However  there  is  a  small  fraction  of  the  photon  signal  that  is  created  due  to 
carrier  diffusion  or  drift.  If  conventional  CL  is  used,  the  actual  distribution  of 
luminescence  is  lost  [2].  In  direct  transport  imaging,  the  light  created  by  the 
recombination  retains  the  spatial  information  of  its  origin.  By  maintaining  the 
actual  distribution  of  the  luminescence,  the  motion  of  the  minority  charge  due  to 
carrier  diffusion  or  drift  can  be  observed. 

B.  SCANNING  ELECTRON  MICROSCOPE  (SEM) 

1.  SEM  Components 

The  two  major  components  of  an  SEM  are  the  electron  column  and  the 
control  console  (Figures  8,  9).  The  electron  column  consists  of  the  electron  gun 
and  various  lenses  that  control  the  beam  path.  The  electron  beam  motion  occurs 
in  an  evacuated  tube  whose  base  pressure  is  maintained  at  ~  - 10'^  Pa  [6]. 


Figure  8.  Electron  Column. 
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Figure  9.  Control  Console. 


The  electron  gun  produces  an  electron  beam  by  thermionic  emission  after 
heating  a  tungsten  filament  to  a  temperature  of  about  2000-2700  K.  It 
accelerates  electrons  in  the  range  0.1-40  keV.  Spot  size  at  the  specimen  is  less 
than  10  nm  at  the  lowest  probe  current  and  penetration  depth  is  ~  about  1-10 
pm  depending  on  the  material  and  the  energy  of  the  incident  beam.  The 
deflection  system  controls  the  magnification.  The  magnification  M  of  the  image  is 
the  ratio  of  the  length  of  the  raster  on  the  viewing  screen  to  the  corresponding 
length  of  the  raster  on  the  specimen  L  (Figure  10).  So  magnification  1000X 
corresponds  to  a  100  pm-wide  raster  on  the  specimen  displayed  on  a  10  cm- 
wide  viewing  screen  [6]. 
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Figure  10.  Deflection  system  inside  the  SEM.  {From  [6]}. 

2.  SEM  Basic  Parameters 

The  electron  probe  size  dp  is  the  diameter  of  the  final  beam  at  the  surface 
of  the  specimen.  The  electron  probe  current  ip  is  the  current  that  impinges  upon 
the  specimen  and  generates  the  various  imaging  signals.  The  electron  probe 
convergence  angle  Op  is  the  half-angle  of  the  cone  of  electrons  converging  onto 
the  specimen.  Finally  the  accelerating  voltage  (Vo)  of  the  electron  beam  defines 
the  electron  beam’s  energy  [6].  All  these  parameters  are  illustrated  in  Figure  11. 


Figure  11.  Major  electron  beam  parameters.  (From  [6]}. 

3.  Electron  Beam-Specimen  Interactions 

If  a  20  kV  accelerating  potential  is  used,  electrons  travel  with  energy  equal 
to  20  keV  with  a  deviation  of  0.5  eV.  Due  to  the  fact  that  the  SEM  column  is 
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under  vacuum  (10-5  Pa),  there  are  very  few  gas  molecules  that  can  collide  with 
electrons  and  deflect  them  out  of  the  beam.  So  almost  all  the  electrons  created 
will  strike  the  surface  of  the  specimen  in  a  circle  of  minimum  diameter  of »  1  nm. 

As  the  beam  electrons  enter  the  specimen,  they  interact  as  negatively 
charged  particles  with  the  specimen  atoms  and  their  corresponding  electrical 
fields.  The  initial  interaction  of  the  electron  beam  with  the  atoms  of  the  specimen 
results  in  an  elastic  scattering  of  the  electrons  with  no  kinetic  energy  loss. 
Electrons,  due  to  this  elastic  scattering,  spread  out  from  the  incident  beam 
footprint.  After  numerous  elastic  scattering  events,  a  portion  of  electrons  finally 
leave  the  specimen  by  a  process  called  “backscattering.”  Simultaneously  with 
elastic  scattering,  beam  electrons  gradually  lose  energy,  transfer  their  energy  to 
specimen  atoms  (inelastic  scattering)  and  secondary  electrons  are  created. 
Backscatter  and  secondary  electrons  are  the  basis  of  the  SEM  imaging  process. 
Additionally,  the  initial  probe  diameter  of  1  nm  has  no  relevance  to  the  image 
resolution.  In  fact,  a  20  keV,  1  nm  diameter  initial  probe,  has  103  times  greater 
linear  dimensions  on  the  entrance  surface  and  109  times  greater  volume 
dimensions  than  the  atomic  layer  under  the  beam  footprint  [6]. 


Figure  12.  Interaction  Volume  for  a  20  keV  beam  striking  Si.  {From  [7]}. 
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In  Figure  12,  the  electron  trajectories  have  been  evaluated  with  the  Monte- 
Carlo  simulation  method.  The  beam  electron  trajectories  that  emerge  as 
backscattered  electrons  are  shown  as  thick  traces. 


All  the  above  can  be  expressed  mathematically  with  the  following 
equations  [6]: 


1.62xl0^"°x 

X  cot^ 

1^7 

12  j 

ds 


keV 
y  cm  J 


=  -7.85x10^"^  In 
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(22) 

(23) 


In  equation  (22)  Q  is  the  probability  of  elastic  scattering  for  angles  greater  than 
the  specified  angle  (po,  Z  is  the  atomic  number  of  the  specimen’s  material  and 

dE 

E  is  the  electron  energy.  In  equation  (23)  —  is  the  rate  of  energy  loss  dE  with 

ds 

distance  traveled  ds ,  p  \s  the  density  of  the  material  in  ,  Ei  is  the  electron 

cm 

energy  at  any  point  in  the  specimen  in  (keV),  A  is  the  atomic  weight  in  — ^ — and 

mole 

J  is  the  average  loss  in  energy  per  event.  The  negative  sign  represents  energy 
loss. 


As  the  beam  energy  increases,  from  equation  (22)  it  can  be  seen  that  the 
electrons  penetrate  more  deeply  into  the  specimen  since  the  probability  of  elastic 
scattering  is  decreased.  From  equation  (23)  it  can  be  seen  also  that  as  energy 
increases,  the  penetration  depth  will  be  larger  since  electrons  enter  the  specimen 
with  more  energy  and  the  rate  of  energy  loss  is  smaller. 

In  order  to  avoid  the  complexity  of  the  3-dimensional  interaction  volume, 
the  latter  can  be  estimated  by  only  one  parameter  which  is  electron  range  and  it 
is  given  by  (due  to  Kanaya  and  Okayama)  [8]; 

0.0276^4  igy 

RKo{pm)  =  ^^ — Eg-  (24) 

Z  p 
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where  A  is  the  atomic  weight  in  — ^ — ,  Z  is  the  atomic  number  of  the 

mole 

specimen’s  material,  p  is  the  density  of  the  material  in  ,  £'o  is  the  beam 

cm 

energy  in  keV  and  Rm  is  the  electron  range  in  pm.  In  this  thesis  the  basic 
semiconductor  material  used  is  GaAs,  so  by  using  equation  (24)  for  evaluating 
the  electron  range  for  a  GaAs  sample  under  the  influence  of  a  beam  with  energy 

Eo  =  3QkeV\ 


Rmipm)  = 


0.0276Z 

7-0.89  ^ 
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The  result  of  equation  (25)  is  valid  if  it  is  assumed  that  the  specimen  is  flat,  thick 
enough  to  be  electron-opaque  and  of  sufficiently  large  lateral  extent  that  there 
are  no  edges  or  boundaries  within  Rko  of  the  beam,  and  the  incident  electron 
beam  is  placed  normal  to  the  specimen  surface.  Due  to  various  uncertainties  of 
the  method,  the  result  should  not  be  stated  beyond  two  significant  figures  [6]. 

4.  SEM  Imaging  Process 

The  information  flow  from  the  SEM  consists  of  the  scan  location  in  x-y 
space  and  a  corresponding  set  of  intensities  from  the  set  of  detectors 
(backscattered  electron,  secondary  electron,  transmitted  electron,  specimen 
current,  cathodoluminescence)  that  monitor  the  beam-specimen  interactions 
simultaneously.  In  this  thesis,  secondary  electron  and  CL  signals  are  collected 
simultaneously. 

a.  Picture  Mode 

In  picture  mode,  the  electron  beam  is  scanned  on  the  specimen 
and  the  resulting  luminescence  pattern  is  recorded  by  the  Charge  Coupled 
Device  (CCD)  camera.  An  example  of  a  picture  of  a  specimen  which  is  going  to 
be  described  later  in  this  thesis  (AIGaAs/GaAs/AIGaAs  heterostructure),  taken  by 
this  mode  of  SEM  operation,  is  shown  in  Figure  13. 
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Figure  13.  Picture  Mode  image  with  application  of  15V  AC  in  a  p-type 
AIGaAs/GaAs/AIGaAs  heterostructure  (874  pm  X  589  pm). 

In  Figure  13,  it  is  interesting  to  observe  the  change  in  the  direction  of  the  electron 
drift  which  is  towards  the  plate  when  the  applied  electric  field  is  positive  (+15V) 
and  outwards  in  the  opposite  case  [negative  electric  field,  (-15V)]. 

b.  Line  Mode 

In  the  line  scan  mode,  the  beam  is  scanned  along  a  single  vector 
on  the  specimen.  With  this  method,  dynamic  electric  fields  can  be  visualized 
better  when  applied  to  relatively  large  areas.  An  example  of  an  image  taken  with 
line  scan  mode  is  shown  in  Figure  14.  In  Figure  14,  the  direction  of  the  electron 
drift  indicates  that  the  plate  (shown  in  red  color  for  identification  purposes  only)  is 
positively  charged.  Also  it  should  be  noted  that  the  larger  spot  in  the  left  side  of 
the  line  scan  is  due  to  the  fact  that  the  SEM  always  holds  the  beam  longer  at  one 
edge  of  the  scanning  area. 


Figure  14.  Line  Scan  Mode  image  with  application  of  15V  DC  in  a  p-type 
AIGaAs/GaAs/AIGaAs  heterostructure  (872  pm  X  293  pm). 
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c.  Spot  Mode 

The  majority  of  the  SEM  imaging,  in  order  to  make  measurements 
of  local  minority  charge  carrier  transport,  has  been  done  with  the  spot  mode  of 
operation.  In  this  mode,  the  electron  beam  is  held  fixed  over  a  specific  point  of 
the  specimen.  An  example  of  this  mode  is  shown  in  Figure  15. 


Figure  15.  Spot  Mode  image  with  application  of  15V  DC  in  a  specific  point 
of  a  p-type  AIGaAs/GaAs/AIGaAs  heterostructure  (216  pm  X 

190  pm). 

The  spot  shown  in  Figure  15  is  the  luminescence  created  by  the  recombination  of 
e-h  pairs  in  response  to  the  interaction  of  the  electron  beam  with  the  specimen. 
Explanation  of  the  shape  of  the  drift  tail  will  be  done  in  the  next  chapter.  All  the 
light  collected  retains  the  spatial  information  of  its  point  of  origin.  The  optical 
image  is  aquired  from  the  CCD  camera  using  “MicroCCD  -  Version  3.13  - 
Diffraction  Limited  -  1999”,  an  image  analysis  software  package.  One  effective 
way  to  study  the  material’s  minority  carrier  properties  is  to  extract  each  “TIFF” 
image  created  in  MicroCCD  to  a  MATLAB  file  and  analyze  the  results  of  the 
intensity  of  luminescence  with  a  method  explained  later  in  this  thesis. 

5.  Laboratory’s  SEM 

In  the  current  work,  the  SEM  used  is  a  JEOL  840A  SEM  with  electron 
energy  range  up  to  40  keV.  In  order  to  observe  luminescence  for  sample 
temperatures  from  5-300  K  there  is  a  liquid  helium-cooled  SEM  stage  (Oxford 

Instruments).  The  optical  detector  is  a  thermoelectrically  cooled  Si  CCD  camera 
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with  a  resolution  of  0.4  |jm/pixel.  The  equipment  is  shown  in  Figure  16  and  an 
overall  schematic  representation  of  the  system  is  given  in  Figure  17. 


MicroCCD 


SEM  electron  Column 


CCD  Camera 


Control  Console 


Figure  16.  JEOL840ASEM. 


Figure  17.  SEM  Overaii  Schematic  Operation  and  Components  {From  [9]}. 


C.  AiGaAs/GaAs/AiGaAs  HETEROSTRUCTURE  SAMPLE 

Figure  18  shows  the  solid  source  molecular  beam  epitaxy  (MBE)-grown 
Alo.9oGao.ioAs/GaAs/Alo.9oGao.ioAs  double  heterostructure,  which  will  be  referred 
to  as  “Sample  #9”  (designed  by  Dr.  Thomas  Boone,  Hitachi  Global  Storage). 
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Figure  18.  AIGaAs/GaAs/AIGaAs  heterostructure  sample. 


The  top  and  bottom  barrier  layers  of  Alo.9oGao.ioAs  were  incorporated  to  reduce 
surface  recombination  effects  and  confine  the  photogenerated  electrons  within 
the  p-type  GaAs  layer.  The  AIGaAs  layers  were  heavily  doped  with  Be 
(Beryllium)  to  Na=2  X  10^®  cm'^  (acceptors/cm^),  with  a  resulted  p-doping  of  the 
GaAs  layer  to  the  mid  1  X  10^®  cm'^.  This  is  known  as  a  modulation  doped 
structure,  where  the  dopant  atoms  are  in  the  AIGaAs  layers  while  the  free  holes 
are  in  the  GaAs,  unaffected  by  impurity  scattering.  A  capping  layer  has  been 
incorporated  on  top  of  the  structure  in  order  to  reduce  surface  oxidization. 
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Additionally,  Van  der  Pauw  mesa  structures  were  etched  into  the  wafer  and  Ti;Au 
nonalloyed  electrical  contacts  were  deposited  and  lithographically  defined  on  top 
of  the  mesas  [10]. 

The  minority  carrier  lifetime  of  the  active  layer  was  measured  with  time 
resolved  photoluminescence  to  be  870  ps.  A  technique  for  imaging  the 
photoluminescent  (PL)  emission  from  the  surface  of  the  semiconductor  under 
increasing  applied  voltage  was  used  by  Dr.  Thomas  Boone  et.al  [10].  The  PL 
displacement  noticed  was  due  to  electron  drift  prior  to  recombination  of  electron- 
hole  pairs.  This  displacement  should  be  (on  average)  equal  to  the  distance  that 
the  electrons  drift  from  their  generation  position  in  their  recombination  lifetime  r  . 
Based  on  the  above  this  distance  is  given  by: 

d  =  VdeT  ^  d  =  {jUeE)r  (26) 

where  d  is  the  distance  traveled  by  electron  while  drifting  from  its  generation 
point,  Vde  is  electron  drift  velocity,  t  is  the  recombination  lifetime,  jue  is  electron 
drift  mobility  and  E  is  the  electric  field  intensity.  Additionally,  for  the  above 

described  p-type  epitaxial  structure,  the  values  used  in  equation  (26)  for 

2 

CtYl 

determining  the  distance  of  electron  drift  were  r«870*10  '^5,  //e«2000 - and 

Vs 

V 

£  max  =  2000 —  (^max  is  the  maximum  value  of  the  electric  field  intensity  used 

cm 

during  the  technique  of  imaging  the  PL  emission  described  above).  Solving  (26) 
with  these  values,  the  resulted  distance  is  d  =  l>5jum  [10]. 

In  Figure  19,  a  part  of  Sample  #9  (top  view),  can  be  seen.  The  red  arrows 
represent  the  places  were  wiring  of  the  sample  has  been  done. 
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Figure  19.  Part  of  Sample  #9  Structure  (3482  pm  X  2611  pm). 

D.  SIMULATIONS 

1.  DC  Modeling 

Fundamental  transport  behavior  can  be  determined  by  imaging  the  charge 
motion  as  a  function  of  applied  electric  field.  Transport  behavior  can  be  predicted 
by  modeling  the  drift  and  diffusion  in  two  dimensions  with  respect  to  generation 
at  a  point  source.  If  an  electric  field  E  is  applied  in  the  x  direction,  the  differential 
equation  that  represents  the  corresponding  intensity  in  the  x  and  y  directions  is 
[2]: 

StU X  L  tU XX  L  tU yy  —  Gt  U7  (27) 

where:  S  =  drift  length,  S  =  jUtE, 

L  =  diffusion  length,  L  =  4^ , 

mx,Gjxx  =  the  first  and  second  derivative  in  x  of  the  carrier  distribution  m , 
vjyy  =  the  second  derivative  in  y  of  the  carrier  distribution  m , 

G  =  generation  rate  of  electron-hole  pairs  [#/(cm^  s)], 

T  =  minority  carrier  (electron)  lifetime. 


31 


The  solution  to  the  differential  equation  (27),  for  a  Gaussian  distribution  e 
of  the  original  beam,  is: 


Intensity  =  — y~2 
27t  L 


□U 

n 


Ko 


s/s^+4L^ 

2Z^ 


sl(x-^y  +{y-r/f 


S(x-?) 

21} 


d^drj 


(28) 


In  equation  (28):  Kq  =  zeroth  order  modified  Bessel  function  of  the  second  kind, 
n  =  parameter  that  defines  the  spread  of  the  generation  region. 

By  simulating  the  diffusion/  drift  behavior  for  the  values  shown  in  Table  3 
below,  the  corresponding  graph  of  Figure  20  represents  the  2D  view  through  the 
distribution  in  x  and  y  directions.  The  color  scale  in  Figure  20  represents  the 
difference  in  intensity  values  derived  from  equation  (28).  The  direction  of  the 
diffusion -drift  is  due  to  the  application  of  an  electric  field  in  the  x-direction  of 
209.50E-4  V/pm.  Figure  20  illustrates  what  is  expected  to  be  seen  in  an  SEM 
image  due  to  the  application  of  a  DC  bias. 


g 

(#/cm^s) 

L 

(pm) 

X 

(range  in  pm) 

y 

(range  in  pm) 

n 

kT 

(eV) 

1 

3.4 

-30:0.5:30 

-30:0.5:30 

0.4 

0.025 

E 

(V/pm) 

PT 

(pm^A/) 

S 

(pm) 

209.50E-4 

pT=eL2/kT 

S=ptE 

Table  3.  Parameters  for  the  evaluation  of  DC  model  for  the  simulation 
of  diffusion/  drift  behavior  in  x  and  y  direction. 


32 


Intensity  Distribution  in  x,  y  Direction 
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Figure  20.  Intensity  distribution  in  x,  y  direction  based  on  simuiation  for 

DC  applied  fieid. 


By  simulating  again  the  diffusion/  drift  behavior  but  now  for  various  values 
of  applied  electric  field  E  and  specifically  for  £'=(0,  30*1 0'"^,  44*10"'^,  70*1 0'"^, 
100*10■^  131. 50*1 0■^  170*10■^  209.50*1 0■^  250*1 0■^  31 3.50*1 0■^  370*1 0■^ 

412*10■^  450*1 0■^  521*10■^  560*1 0■^  600*10'^)  V/pm,  the  corresponding  graph 
is  shown  in  Figure  21  below.  Additionally,  the  other  parameters  needed  for  the 
evaluation  of  (29)  were  defined  as  shown  in  Table  4. 
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g 

(#/cm^s) 

L 

(pm) 

X 

(range  in  pm) 

y 

(range  in  pm) 

n 

kT 

(eV) 

1 

3.4 

-30:0.5:30 

0 

0.4 

0.025 

E 

PT 

S 

(V/pm) 

(pm^A/) 

(pm) 

values  mentioned  above 

pT=eL^/kT 

S=ptE 

Table  4.  Parameters  for  the  evaluation  of  DC  model  for  the  simulation 

of  diffusion/  drift  behavior. 
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Figure  21.  Simulation  for  varying  values  of  DC  applied  fieid. 


The  results  shown  in  Figure  21  reflect  the  increasing  drift  length  and 
transport  of  charge  with  increasing  electric  field.  It  can  be  seen  also  that,  as  the 
applied  voltage  increases,  the  peak  intensity  of  each  curve  decreases  and  shifts 
to  the  right  (with  respect  to  the  central  value  x=0).  All  the  areas  under  the  curves 
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are  equal.  This  can  be  shown  by  use  of  the  trapezoid  rule  for  evaluating  the  area 
under  each  curve  (numerical  analysis).  If  the  above  model  is  simulated  for  a 
much  greater  number  of  electric  field  values  and  the  value  of  peak  intensity  is 
defined  for  each  case  (different  applied  voltage),  a  plot  of  normalized  intensity 
(with  respect  to  E=0  V/pm)  versus  applied  electric  field  can  be  created  (Figure 
22).  In  Figure  23,  a  regression  analysis  has  been  made,  in  order  to  define  the 
curve  fitting  parameters  of  the  line  connecting  the  normalized  peak  intensity 
values  with  respect  to  applied  electric  field.  It  should  be  noted  that  the  axes  in 
Figure  23  have  been  inverted  only  for  simplifying  the  fitting  parameter  procedure. 
The  MATLAB  code  used  for  the  generation  of  the  current  graphs  can  be  found  in 
Appendix  A. 


Figure  22.  Intensity  ratio  vs.  applied  electric  field  based  on  simulation  for 

DC  applied  field. 
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2.  AC  Modeling 

Similar  to  the  previous  technique  of  obtaining  the  DC  field  behavior, 
an  analogous  method  for  the  AC  field  has  been  developed.  The  model  for  the  AC 
field  is  the  following: 


Intensity  =  |  ^  |  j 


^JSitf+4L^ 

2L^ 


S{t){x-S) 
21} 


d^drjdt 


where:  S  =  drift  length,  S{t)  =  i^xE{t) ,  (29) 

T  =  diffusion  length,  L  =  4^ , 

G  =  generation  rate  of  electron-hole  pairs  [#/(cm3  s)], 

T  =  minority  carrier  (electron)  lifetime, 

E{t)  =  Asm{cot)  =  sinusoidal  AC  electric  field, 

a)  =  lnf  where  /  is  the  frequency  used  for  the  AC  current.  The  rest  of 
parameters  are  the  same  as  for  equation  (28). 
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By  simulating  the  diffusion/  drift  behavior  for  the  values  shown  in  Table  5 
below,  the  corresponding  graph  of  Figure  24  represents  the  2D  view  of  the 
distribution  in  x  and  y  directions  for  the  case  of  AC  bias  (for  t=2s).  The  color  scale 
in  Figure  24  represents  the  difference  in  intensity  values  derived  from  equation 
(29).  The  distribution  is  symmetric  due  to  the  application  of  AC  bias. 
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(range  in  pm) 
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(range  in  pm) 

n 

kT 

(eV) 
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E(t) 

A 

PT 

S(t) 

f 

to 

(V/pm) 

(V/pm) 

(pm^A/) 

(pm) 

(Hz) 

(rad/s) 

E(t)=Asin(u)t) 

209.50E-4 

pT=eL2/kT 

S(t)=pTE 

333 

U)=2TTf 

Table  5.  Parameters  for  the  evaluation  of  AC  model  for  the  simulation 
of  diffusion/  drift  behavior  in  x  and  y  direction. 

Intensity  Distribution  in  x,y  Direction 
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Figure  24.  Intensity  distribution  in  x,  y  direction  based  on  simulation  for 

AC  applied  field. 

37 


By  simulating  the  diffusion/  drift  behavior  for  various  values  of  applied 
electric  field  E  and  specifically  for  amplitude  of  electric  field  (peak-to-peak), 
A=(0,  44*1 0■^  131 .50*1 0■^  209.50*1 0■^  300*1 0■^  412*10'^)  V/pm,  the 
corresponding  graph  is  shown  in  Figure  25.  Additionally,  the  other  parameters 
needed  for  the  evaluation  of  (30)  were  defined  as  shown  in  Table  6. 
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(range  in  pm) 
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(range  in  pm) 
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kT 

(eV) 
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PT 
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f 

U) 
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as  above 
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Table  6.  Parameters  for  the  evaluation  of  AC  model  for  the  simulation 

of  diffusion/  drift  behavior. 
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Figure  25.  Simulation  for  varying  vaiues  of  AC  appiied  fieid. 

38 


The  results  shown  in  Figure  25  reflect  the  increasing  drift  length  and 
transport  of  charge  with  increasing  electric  field.  The  value  of  the  peak  intensity 
of  each  curve  is  directly  dependent  on  time  (by  increasing  the  time  the  value  of 
the  intensity  increases  too)  but  by  evaluating  the  ratio  l/lo,  the  results  (normalized 
intensity  ratios)  are  independent  of  time.  This  can  be  seen  clearly  in  Figure  26 
where  a  comparison  of  normalized  intensity  ratios  with  respect  to  different  times 
(2  sec  and  5  sec)  has  been  performed.  In  Figure  25,  it  can  also  be  seen  that,  as 
the  applied  voltage  increases,  the  peak  intensity  of  each  curve  decreases  but 
remains  as  expected,  in  the  center  of  the  graph  (with  respect  to  the  central  value 
x=0).  All  the  areas  under  the  curves  are  equal.  This  can  be  shown  again  by  use 
of  the  trapezoid  rule  for  evaluating  the  area  under  each  curve  (numerical 
analysis).  If  the  values  of  peak  intensity  are  defined  for  each  case  (different 
applied  voltage),  a  plot  of  normalized  intensity  (with  respect  to  E=0  V/pm)  versus 
applied  electric  field  is  shown  in  Figure  27  below.  Also  in  Figure  28,  a  regression 
analysis  has  been  made,  in  order  to  define  the  curve  fitting  parameters  of  the  line 
connecting  the  normalized  peak  intensity  values  with  respect  to  applied  electric 
field.  Note  that  the  axes  in  Figure  28  have  been  inverted  only  for  simplifying  the 
fitting  parameter  procedure. 


Figure  26.  Comparison  of  normalized  intensity  ratios  for  different  times. 
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Figure  27. 


Normalized  Peak  Intensity  Ratio  (Theory)  for  AC  Model  (t=5s) 


E-Field  (V/nm) 


Intensity  ratio  vs.  applied  electric  field  based  on  simulation  for 

AC  applied  field. 


It  should  be  noted  at  this  point  that  the  values  of  the  electric  field  shown  in 
Figures  27  and  28  are  the  root-mean -square  values  (rms  or  effective  values)  that 
are  derived  by  (for  a  pure  sinusoidal  waveform): 
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/  (x)  =  sin  X 


rms  =  foA  =  f^{x)dx  =  j  sin^  xdx  =  =  A =  0.707  x  A 
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sin  An  sin  0 

- 1 - 

4  4 


=  n 


(30) 


where  A  is  the  peak-to-peak  amplitude  of  the  AC  bias. 

Finally  in  Figure  29,  a  comparison  of  the  normalized  intensity  ratios  based 
on  the  DC  and  AC  simulations  for  a  bias  of  209.50E-4  V/pm,  is  illustrated.  The 
peak  intensity  with  AC  bias  is  higher  than  the  equivalent  with  DC  bias,  reflecting 
a  slightly  lower  effective  field. 


Figure  29.  Comparison  of  DC  and  AC  simulation  models  for  E-field 

amplitude  of  209.50E-4  V/pm. 
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IV.  ELECTRIC  FIELD  DISTRIBUTIONS:  SIMULATION  AND 

MEASUREMENT 

A.  2-D  ELECTRIC  FIELD  SIMULATIONS 

In  order  to  determine  the  extent  to  which  the  charge  imaging  can  be  used 
to  map  electric  fields  with  high  resolution,  charge  motion  is  imaged  as  a  function 
of  applied  electric  field  in  a  variety  of  planar  contact  geometries.  Transport 
behavior  under  the  combined  influence  of  both  diffusion  and  drift  can  be 
predicted  by  using  the  simulation  models  described  in  Chapter  III  for  diffusion 
and  drift  in  2D  following  generation  at  a  point  source.  The  overall  method 
described  in  this  chapter  is  based  on  a  comparison  of  the  experimental  field 
measurements  to  2D  finite  element  field  simulations  for  the  structures  of  interest. 

All  the  following  cases  have  been  evaluated  experimentally  on  the 
Alo.9oGao.ioAs/GaAs/Alo.9oGao.ioAs  double  heterostructure  (Sample  #9)  described 
in  detail  in  the  previous  chapter. 

1.  Case  1 

In  Figure  30  an  overview  of  Sample  #9  is  shown  to  illustrate  the  overall 
procedure.  In  Case  1,  15  V  DC  has  been  applied  across  the  contacts  as  shown 
in  Figure  30.  The  SEM  is  operated  in  spot  mode.  The  location  of  the  electron 
beam  with  respect  to  the  sample  is  represented  with  a  red  circle  (for  each 
position,  PI,  P2,  Pn).  The  probe  current  was  3X10-10  A  and  the  electron  beam 
energy  has  been  fixed  at  30  keV.  The  dimensions  of  the  study  area  also  can  be 
seen  in  Figure  30. 
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Figure  30.  Sample  #9  -  Case  1. 


In  order  to  normalize  the  data  from  each  measurement,  two  sets  of  data 
are  obtained  in  each  position:  one  with  0  V  and  one  with  the  applied  15V  DC. 
The  procedure  of  locally  mapping  the  electric  field  based  on  the  direct  transport 
imaging  technique  includes  the  following: 

•  Image  the  charge  recombination  length  in  each  position  first  with  no 
bias  and  then  with  applied  bias. 

•  Extract  the  “.tif  images  created  in  the  MicroCCD  software  to  a 
MATLAB  program  and  find  the  value  of  peak  intensity  in  each  case. 

•  Evaluate  the  normalized  intensity  of  each  case  based  on  the 
measurements  with  and  without  bias. 

•  Create  a  graph  of  normalized  intensity  ratio  as  a  function  of 
distance,  based  on  the  experimental  measurements. 

•  Combine  the  previous  graph  with  the  data  obtained  from  the 
theoretical  model  of  intensity  as  a  function  of  distance  for  DC 
applied  electric  field  as  analyzed  previously  and  apply  the  curve 
fitting  equation  to  obtain  a  graph  of  electric  field  as  a  function  of 
distance. 

The  MATLAB  code  created  for  this  purpose  can  be  seen  in  Appendix  B. 
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The  results  of  the  above  procedure  for  DC  applied  field  can  be  seen  in 
Figures  31  and  32.  Figure  31  represents  the  normalized  intensity  ratio  with 
respect  to  the  distance,  based  on  the  SEM  experimental  measurements.  Figure 
32  represents  the  electric  field  magnitude  with  respect  to  distance,  based  on  the 
procedure  developed  in  order  to  map  the  electric  field  through  the  transport 
imaging  technique.  By  comparing  the  two  figures  it  can  be  seen  that  when  the 
normalized  intensity  ratio  increases,  the  magnitude  of  the  electric  field  at  the 
same  point  decreases  and  vice-versa. 


Figure  31.  Normalized  intensity  ratio  vs.  Distance  for  Case  1  (DC  applied 

field). 


E-Field  vs  Distance 


Figure  32.  E-Field  vs.  Distance  for  Case  1  (DC  applied  field). 
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By  repeating  the  above  method  with  AC  applied  electric  field  (15  V  AC) 
and  keeping  all  the  other  parameters  unchanged,  the  corresponding  graphs  can 
be  seen  in  Figures  33  and  34.  The  reason  for  using  AC  electric  field  is  for 
comparison  purposes  with  respect  to  DC  applied  field  and  for  increasing  image 
information  due  to  the  fact  that  the  drift  and  diffusion  is  simultaneously  observed 
in  both  directions. 


Experimental  Measurement 
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Figure  33.  Normalized  intensity  ratio  vs.  Distance  for  Case  1  (AC  applied 

field). 


E-Field  vs  Distance 


Distance  (mm) 


Figure  34.  E-Field  vs.  Distance  for  Case  1  (AC  applied  field). 
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From  Figures  33  and  34,  it  can  be  seen  that  the  application  of  AC  bias 
instead  of  DC  bias  does  not  alter  either  the  shape  of  the  graphs  or  significantly 
the  magnitude  of  the  electric  field  at  the  same  points  of  interest.  This  can  be  seen 
clearly  in  Figure  35  where  a  comparison  of  the  AC-DC  applied  field  results  has 
been  done.  In  the  DC  case  the  magnitude  of  the  electric  field  is  higher  on  the 
regions  where  the  field  is  stronger  and  lower  than  the  AC  values  in  regions  where 
the  field  is  lower. 
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Figure  35.  Comparison  of  the  results  of  magnitude  of  electric  field  with 

DC  and  AC  bias  for  Case  1. 
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2.  Case  2 

In  Case  2,  15  V  DC  has  been  applied  between  the  contacts  as  shown  in 
Figure  36.  The  SEM  is  again  operated  in  spot  mode.  The  location  of  the  electron 
beam  with  respect  to  the  sample  is  represented  with  a  red  circle  (for  each 
position).  The  probe  current  was  3X10'^°  A  and  the  electron  beam  energy  was  30 
keV. 


Note  that  Figure  36  demonstrates  the  electric  potential  distribution  as 
calculated  in  COMSOL  MULTIPHYSICS  (Version  3.3)  software.  The  magnitude 
of  the  potential  in  each  region  is  represented  by  different  colors  and  the 
corresponding  values  are  shown  in  the  color  scale  in  the  right  part  of  the  graph. 
With  application  of  15  V  DC,  the  corresponding  normalized  intensity  ratio  and 
electric  field  graphs  are  shown  in  Figures  37  and  38. 
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Exp>erimental  Measurement 


Figure  37.  Normalized  intensity  ratio  vs.  Distance  for  Case  2  (DC  applied 

field). 
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Figure  38.  E-Field  vs.  Distance  for  Case  2  (DC  applied  field). 
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For  the  AC  case,  the  corresponding  normalized  intensity  ratio  and  electric 
field  graphs  are  shown  in  Figures  39  and  40. 


Figure  39.  Normalized  intensity  ratio  vs.  Distance  for  Case  2  (AC  applied 

field). 
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Figure  40.  E-Field  vs.  Distance  for  Case  2  (AC  applied  field). 
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As  in  Case  1,  the  comparison  of  DC  and  AC  applied  bias  results  show 
(Figure  41)  that  the  DC  values  are  higher  than  the  AC  results  (almost  20%  in  the 
strong  field  region)  and  a  little  lower  than  the  AC  in  the  weak  electric  field  region 
(almost  3%). 
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Figure  41 .  Comparison  of  the  results  of  magnitude  of  electric  field 

DC  and  AC  bias  for  Case  2. 
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3.  Case  3 

In  Case  3,  15  V  DC  has  been  applied  between  the  contacts  as  shown  in 
Figure  42  with  the  same  experimental  parameters  as  in  the  previous  cases.  For 
the  case  of  DC  bias,  the  corresponding  normalized  intensity  ratio  and  electric 
field  graphs  are  shown  in  Figures  43  and  44. 


Figure  42.  Sample  #9  -  Case  3. 
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Figure  43.  Normalized  intensity  ratio  vs.  Distance  for  Case  3  (DC  appiied 

field). 
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E-Field  vs  Distance 


Distance  (mm) 

Figure  44.  E-Field  vs.  Distance  for  Case  3  (DC  applied  field). 

Additionally,  a  similar  experimental  measurement  was  made  with  the 
same  parameters,  except  that  the  cabling  and  the  polarity  have  been  reversed. 
The  result  is  illustrated  in  Figure  45  below.  The  corresponding  normalized 
intensity  ratio  and  electric  field  graphs  are  shown  in  Figures  46  and  47  below. 


Figure  45.  Sample  #9  (Reversed  cabling  and  polarity). 
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Figure  46.  Normalized  intensity  ratio  vs.  Distance  for  Case  3  (DC  applied 

field,  reversed  cabling  and  polarity). 
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Figure  47.  E-Field  vs.  Distance  for  Case  3  (DC  applied  field,  reversed 

cabling  and  polarity). 
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Although  the  cases  shown  in  Figures  42  and  45  might  seem  different,  they 
are  identical.  In  both  cases  the  electric  potential  is  the  same  as  the  polarity  of  the 
applied  bias  and  the  SEM  cabling  to  the  sample  have  been  reversed 
simultaneously.  By  comparing  Figures  44  and  47  it  can  be  seen  that  due  to  the 
reversed  cabling  and  polarity  in  the  same  structure,  the  magnitude  of  the  electric 
field  with  respect  to  distance  at  each  point  of  interest,  is  not  identical  as 
theoretically  was  anticipated.  The  reason  for  the  same  electric  potential  for  both 
cases  is  due  to  the  fact  that  the  reversed  polarity  (from  positive  +15V  to  negative 
-15V)  reverses  the  field  but  by  switching  the  cabling  of  the  SEM  the  resulting  field 
is  the  initial  case.  This  can  be  seen  clearly  in  Figure  48  where  a  comparison  of 
the  two  cases  is  demonstrated.  The  voltage  drop  at  metal-semiconductor 
material  interface  is  not  always  symmetric.  For  this  reason  there  are  some 
fluctuations  in  the  magnitude  of  the  electric  field  {~4%). 
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Figure  48.  Comparison  for  Case  3  resuits  (DC  appiied  field,  reversed 

cabling  and  polarity). 
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4.  Case  4 

In  Case  4,  15  V  DC  has  been  applied  between  the  contacts  as  shown  in 
Figures  49  and  50. 


Figure  49.  Sample  #9  (SEM  photo)  -  Case  4  (2031  pm  X  2090  pm). 
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Figure  50.  Sample  #9  -  Case  4  (Created  in  COMSOL  MULTIPHYSICS). 
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The  SEM  is  operated  in  spot  mode.  Twelve  measurements  of  the  intensity  have 
been  made  across  the  yellow  arrow.  The  probe  current  was  3X10'^°  A  and  the 
electron  beam  energy  was  at  30  keV.  The  corresponding  graphs  of  the 
normalized  intensity  ratio  and  electric  field  are  shown  in  Figures  51  and  52. 
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Figure  51.  Normalized  intensity  ratio  vs.  Distance  for  Case  4  (DC  applied 

field). 
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Figure  52.  E-Field  vs.  Distance  for  Case  4  (DC  applied  field). 
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By  applying  15  V  AC  in  the  above  configuration  the  corresponding  graphs 
of  the  normalized  intensity  and  electric  field  are  given  in  Figures  53  and  54. 
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Figure  53.  Normalized  intensity  ratio  vs.  Distance  for  Case  4  (AC  applied 

field). 


Figure  54.  E-Field  vs.  Distance  for  Case  4  (AC  applied  field). 
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From  Figure  55  it  can  be  seen  that  the  DC  results  for  the  magnitude  of  the 
electric  field  are  higher  than  the  equivalent  of  the  AC  results,  as  expected,  by 
about  21%.  Again  in  the  region  of  weak  electric  field  the  AC  and  DC  values 
converge. 


Figure  55.  Comparison  of  the  results  of  magnitude  of  electric  field  with 

DC  and  AC  bias  for  Case  4. 
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5.  Case  5 

In  Case  5,  15  V  DC  have  been  applied  between  the  contacts  as  shown  in 
Figures  56  and  57. 


Figure  56.  Sample  #9  (SEM  photo)  -  Case  5  (2031  pm  X  2090  pm). 
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Sample  #9  -  Case  5  (Created  in  COMSOL  MULTIPHYSICS). 
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The  SEM  is  operated  in  spot  mode.  Eighteen  consecutive  measurements 
of  the  intensity  have  been  made  across  the  yellow  arrow.  The  probe  current  was 
3X10"^°  A  and  the  electron  beam  energy  was  30  keV.  The  corresponding  graphs 
of  the  normalized  intensity  ratio  and  the  electric  field  are  shown  in  Figures  58  and 
59. 
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Figure  58.  Normalized  intensity  ratio  vs.  Distance  for  Case  5  (DC  applied 

field). 


E-Field  vs  Distance 


Figure  59.  E-Field  vs.  Distance  for  Case  5  (DC  applied  field). 
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By  applying  15  V  AC  in  the  above  configuration  the  corresponding  graphs 
of  normalized  intensity  ratio  and  electric  field  are  shown  in  Figures  60  and  61 
below. 
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Figure  60.  Normalized  intensity  ratio  vs.  Distance  for  Case  5  (AC  applied 

field). 
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Figure  61.  E-Field  vs.  Distance  for  Case  5  (AC  applied  field). 
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In  Figure  62  it  can  be  seen  again  that  the  AC  results  for  the  electric  field 
magnitude  are  lower  than  the  equivalent  DC  results  by  about  20%  except  in  the 
regions  of  weak  electric  field  where  the  AC  values  are  almost  identical  with  the 
equivalent  DC  values. 
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Figure  62.  Comparison  of  the  results  of  magnitude  of  electric  field  with 

DC  and  AC  bias  for  Case  5. 
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6.  Case  6 

In  Case  6,  15  V  DC  has  been  applied  between  the  contacts  as  shown  in 


Figures  63  and  64. 


Figure  63.  Sample  #9  (SEM  photo)  -  Case  6  (2031  pm  X  2090  pm). 


Figure  64.  Sample  #9  -  Case  6  (Created  by  COMSOL  MULTIPHYSICS). 
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The  SEM  is  operated  in  spot  mode.  Seventeen  consecutive 
measurements  of  the  intensity  have  been  made  across  the  yellow  arrow.  The 
probe  current  was  3X10'^°  A  and  the  electron  beam  energy  was  at  30  keV.  The 
electric  field  between  the  plates  is  uniform  except  in  the  region  right  next  to  the 
plates.  The  corresponding  graphs  can  be  seen  below  in  Figures  65  and  66. 
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Figure  65.  Normalized  intensity  ratio  vs.  Distance  for  Case  6  (DC  applied 

field). 


E4^ield  vs  Distance 


0.095 


0.09  - 


0.085 


0.08  - 


0.075  - 


0.07 


o 

"  o 

- 1 

1 

1 

1 

1 

— 

^ . 

o 

c 

)  O 

o 

. O"" 

. <: 

o 

o 

o 

o 

o 

'D . 

_ 

_ 

_ 

_ 

_ 

1 

_ 

0.05  0.1  0.15  0.2  0.25  0.3  0.35 

Distance  (mm) 


Figure  66.  E-Field  vs.  Distance  for  Case  6  (DC  applied  field). 
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By  applying  15  V  AC  in  the  above  configuration  the  corresponding  graphs 
are  shown  in  Figures  67  and  68  below. 
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Figure  67.  Normalized  intensity  ratio  vs.  Distance  for  Case  6  (AC  bias). 


E-Field  vs  Distance 


Figure  68.  E-Field  vs.  Distance  for  Case  6  (AC  bias). 


In  Figure  69  below  it  can  be  seen  that  the  values  of  the  electric  field  with 
the  AC  bias  case  are  lower  than  the  equivalent  with  DC  bias  in  the  whole  region 
by  about  25%. 
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The  fluctuations  of  the  magnitude  of  the  electric  field  in  the  region  between 
the  contacts  seem  to  be  different  from  what  was  expected  since  the  field  is 
uniform  except  for  a  small  area  in  proximity  to  the  metal  contacts.  But  these 
fluctuations  are  very  small  in  magnitude  and  this  can  be  seen  clearly  in  Figure  70 
where  the  y-axis  has  been  expanded  to  be  consistent  with  earlier  graphs  and  the 
values  of  the  electric  field  are  almost  all  in  a  line  as  expected  for  a  uniform 
electric  field.  However,  the  experimental  technique  is  sensitive  to  any  variations 
in  field  associated  with  material  non-uniformity. 
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Figure  69.  Comparison  of  the  results  of  magnitude  of  eiectric  field  with 

DC  and  AC  bias  for  Case  6. 
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Figure  70.  DC  and  AC  bias  comparison  (different  scaiing). 
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By  comparing  all  the  previous  cases,  the  consistent  variation  between  the 
DC  and  AC  applied  field  cases  suggests  either: 

•  Some  additional  correction  is  needed  to  the  AC  modeling,  or 

•  There  is  an  increased  contact  resistance  when  AC  bias  is  applied. 
The  fact  that  the  difference  depends  on  the  magnitude  of  the  E-field  suggests  a 
possible  contact  effect,  even  at  this  relatively  low  frequency. 

B.  VISUALIZATION  OF  THE  ELECTRIC  FIELD 

One  very  important  aspect  of  the  method  described  in  this  thesis  is  the 
capability  of  visualizing  the  direction  of  the  electric  field  locally.  In  order  to 
demonstrate  this,  the  case  shown  in  Figure  71  will  be  used.  The  SEM  is  operated 
in  line  mode  (scanning  direction  is  represented  with  the  yellow  arrow).  The 
results  are  illustrated  in  Figure  72. 


Figure  71.  Direction  of  the  electric  field  with  SEM  in  line  mode  (874  pm  X 

589  pm). 
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Figure  72.  Line  mode  in  defining  the  direction  of  the  electric  field  (each 

image  is  872  |jm  X  293  |jm). 

In  Figure  72  the  direction  of  the  electric  field  can  be  directly  observed.  In 
part  (a)  of  Figure  72,  no  voltage  is  applied.  In  part  (b)  15  V  DC  have  been  applied 
between  the  contacts.  The  drift  of  the  minority  carriers  (electrons)  is  towards  the 
positively  charged  plate.  In  part  (c)  15  V  AC  have  been  applied  to  the  above 
configuration.  This  results  in  a  symmetrical  shape  of  the  line  mode  as  the  voltage 
polarity  oscillates.  Again,  it  should  be  noted  that  the  bright  spot  shown  in  the  left 
of  each  line  is  an  artifact  of  the  SEM  imaging  process  due  to  the  synchronization 
of  the  electron  beam  (SEM  always  scans  the  beam  longer  along  one  edge  of  the 
scanning  area). 

Additionally,  by  operating  the  SEM  in  picture  mode,  generation  over  the 
whole  region  is  illustrated  in  Figure  73. 
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Figure  73.  Various  cases  of  direction  of  diffusion  length  according  to  the 
applied  voltage  (each  image  is  874  pm  X  589  pm). 


In  part  (a)  of  the  figure  no  voltage  is  applied.  In  part  (b)  -15  V  DC  bias  has 
been  applied  to  the  contact.  Minority  carriers  drift  opposite  to  the  negatively 
charged  contact.  In  part  (c)  +15  V  DC  has  been  applied.  In  this  case,  minority 
carriers  drift  towards  the  positively  charged  plate.  Finally,  in  part  (d)  of  the  figure, 
15  V  AC  has  been  applied.  This  causes  minority  carriers  to  drift  towards  and 
away  from  the  plate  according  to  the  sign  (+  or  -)  of  the  instantaneous  voltage. 
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Finally,  by  operating  SEM  in  spot  mode,  case  5  (paragraph  A-5),  is  used 
to  visualize  the  E-field  vector  and  to  illustrate  the  differences  between  the  DC 
and  AC  imaging  process. 


Figure  74.  Various  positions  of  electron  beam  (spot  mode)  in  Sample  #9 


(874  pm  X  589  pm). 
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Figure  75.  Direction  of  the  electric  field  in  case  5  based  on  the  positions 


specified  in  Figure  74  for  DC  bias  (each  image  is  266  pm  X  239 

pm). 

In  Figure  75,  images  of  the  charge  motion  for  varying  directions  of  the 
electric  field  are  shown  based  on  the  generation  spots  of  the  electron  beam  as 
specified  in  Figure  74.  Sample  #9  is  a  p-type  semiconductor,  so  the  minority 
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carriers  are  electrons.  The  consecutive  images  in  Figure  75  represent  the  drift 
motion  of  minority  carriers  (electrons)  and  this  can  be  seen  clearly  from  the 
direction  of  the  electric  field  based  on  the  contacts;  electrons  follow  the  electric 
field  lines  created  from  the  potential  difference  between  the  two  contacts. 
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Figure  76.  Direction  of  the  electric  field  in  case  5  based  on  the  positions 
specified  in  Figure  74  for  AC  bias  (each  image  is  266  pm  X  239 

pm). 


In  Figure  76,  the  same  imaging  technique  has  been  adopted  with  the 
difference  that  AC  bias  has  been  applied  as  in  the  second  part  of  case  5.  The  AC 
bias  creates  a  symmetry  in  the  drift  behavior  but  again  the  drift  is  towards  the 
positive  contact  as  the  polarity  oscillates. 

Based  on  the  analysis  that  has  been  done  in  the  current  chapter  it  can  be 
seen  that  the  charge  imaging  technique  is  a  powerful  “tool”  to  map  electric  fields 
with  high  resolution.  By  imaging  the  motion  of  the  minority  carriers  (electrons)  as 
a  function  of  applied  electric  field  in  a  variety  of  planar  contact  geometries,  both 
the  magnitude  and  the  direction  of  the  electric  field  vector  can  be  defined  with 
high  accuracy. 

C.  COMSOL  MULTIPHYSICS 

COMSOL  MULTIPHYSICS  is  a  finite  element  analysis  and  solver  software 
package  for  various  physics  and  engineering  applications,  especially  coupled 
phenomena  or  multiphysics.  COMSOL  MULTIPHYSICS  also  offers  an  extensive 
and  well-managed  interface  to  MATLAB  and  its  toolboxes  for  a  large  variety  of 
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programming,  preprocessing  and  postprocessing  possibilities.  A  similar  interface 
is  offered  to  COMSOL  Script.  The  packages  can  be  used  with  every  operating 
system  (MS  Windows,  Mac,  Linux,  Unix).  In  addition  to  conventional  physics- 
based  user-interfaces,  COMSOL  MULTIPHYSICS  also  allows  for  entering 
coupled  systems  of  partial  differential  equations  [11]. 

For  the  purposes  of  the  current  thesis,  the  Electromagnetics  Module  of 
Version  3.2b  and  the  AC/DC  Module  of  Version  3.3  have  been  used.  All  the 
semiconductor  planar  structures  (portions  of  Sample  #9)  used  have  been 
designed  with  the  CAD  capability  of  the  software.  This  particular  module  solves 
problems  in  the  general  areas  of  electrostatic  fields,  magnetostatic  fields,  and 
quasi-static  fields.  The  equations  for  electromagnetics  are  automatically  available 
in  all  of  the  application  modes,  making  simultaneously  the  nonstandard  modeling 
easily  accessible.  The  electrostatics  (generalized)  part  of  the  AC/DC  Module 
used  simulates  electric  fields  and  currents  in  dielectric  and  conductive  materials 
[11]. 

Finally,  all  the  simulations  of  the  electric  fields  in  COMSOL  have  been 
done  in  2D.  The  modeling  in  2D  usually  represents  some  3D  geometry  under  the 
assumption  that  nothing  changes  in  the  third  dimension.  In  order  to  use  the  2D 
cross-section  view  of  a  problem,  there  are  two  approaches:  the  first  is  when  it  is 
assumed  that  there  is  no  variation  of  the  solution  in  one  particular  dimension. 
The  second  is  when  the  influence  of  the  finite  extension  in  the  third  dimension 
can  be  neglected.  For  this  reason,  in  the  simulations  demonstrated  in  the  next 
chapter,  even  though  the  2D  geometry  is  used,  there  are  calculations  of  the  z- 
component  of  the  electric  field  [1 1]. 
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V.  APPLICATIONS  OF  THE  TRANSPORT  IMAGING 

TECHNIQUE 


A.  CONTACT  RESISTANCE  MEASUREMENTS 

When  a  voltage  is  applied  via  contacts  to  a  semiconductor  structure,  there 
is  some  voltage  drop  associated  with  the  semiconductor-metal  interface.  The 
magnitude  of  this  voltage  reduction  depends  on  the  quality  and  nature  of  the 
contact’s  material.  A  standard  approach  for  determining  the  contact  resistance  is 
a  4-point  measurement.  The  disadvantage  of  the  method  is  the  necessity  to  use 
two  current  carrying  probes  and  two  voltage  measuring  probes. 

Transport  imaging  offers  a  direct  way  to  evaluate  the  voltage  drop  at  the 
contacts  with  a  2-point  measurement.  The  advantage  of  the  method  is  that  no 
additional  processing  is  required.  Since  the  structure  of  Sample  #9  is  based  on 
Schottky  diode  contact  formation,  deviations  from  the  ohmic  behavior  are 
expected.  For  this  reason,  diode  l-V  characteristics  are  going  to  be  explored  in 
the  current  chapter. 

1.  Case  1 

For  this  case  the  configuration  is  shown  in  Figure  77. 


In  Figure  77,  15  V  DC  has  been  applied  to  the  outer  plates  and  the 
distance  between  the  inner  plates  (edge  to  edge)  is  185  pm.  The  electron  beam 
is  represented  with  the  red  circle  (not  to  scale).  More  specifically,  the  area  of 
interest  is  illustrated  in  Figure  78  by  using  the  COMSOL  software. 


Figure  78.  COMSOL  simulation  of  the  electric  field  potential  for  Case  1. 


In  Figure  78,  the  electric  potential  is  represented  with  color  variations 
defined  in  the  corresponding  scale  in  the  right  part  of  the  figure  above.  The 
contour  lines  represent  equipotential  lines.  The  red  circle  shows  the  position  of 
the  electron  beam  (not  to  scale).  The  x,  y  and  z  components  of  the  generated 
electric  field  between  the  metal  contacts  which  are  superjacent  to  the 
semiconductor  material  (Schottky  junction),  can  be  seen  in  the  following  Figures 
79,  80  and  81. 
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Figure  79. 


Figure  80. 


COMSOL  simulation  of  the  electric  field  (x-component)  for 

Case  1. 


ov 

15V 

tm 


COMSOL  simulation  of  the  electric  field  (y-component)  for 

Case  1. 
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Figure  81 .  COMSOL  simulation  of  the  electric  field  (z-component)  for 

Case  1. 


It  should  be  noted  again  that  even  though  the  model  used  in  COMSOL 
was  based  on  a  2D  geometry,  the  z  component  has  the  physical  meaning  that 
the  magnitude  of  the  electric  field  in  the  z-direction  is  constant  so  the  direction  of 
the  field  is  tangent  to  the  x-y  plane.  The  generated  graphs  for  the  electric  field, 
along  the  center  line  between  the  inner  plates,  in  each  direction  can  be  seen  in 
Figures  82,  83  and  84  below.  The  procedure  used  in  order  to  generate  the 
following  graphs  is: 

•  In  COMSOL,  a  perpendicular  line  is  plotted  between  the  inner 
plates,  passing  through  the  position  of  the  electron  beam. 

•  COMSOL  automatically  generates  a  graph. 

•  Extract  the  data  in  a  “.txt”  file. 

•  Import  the  data  into  a  MATLAB  code  (Appendix  C)  and  generate 
the  corresponding  plots. 
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x  10^  Electric  Field  -  x  component 


Figure  82.  Electric  field  (x-component)  for  Case  1. 


Figure  83.  Electric  field  (y-component)  for  Case  1. 
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Figure  84.  Electric  field  (z-component)  for  Case  1. 


In  Figure  85,  the  total  electric  field  from  the  simulation  is  presented.  It  can 
be  seen  that  there  is  a  perfect  symmetry  in  the  graph  and  this  is  exactly  as 
expected  since  the  model  assumes  completely  uniform  material  and  contact 
parameters. 
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Figure  85.  Total  electric  field  for  Case  1. 
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By  determining  the  value  of  the  electric  field  in  the  middle  of  the  distance 
between  the  plates  (i.e.,  185/2=92.5  pm),  in  each  direction,  the  resultant  value  is; 

£.  =  -5.57x10"—  ,  £v  =  -64.8— ,  £z  =  5.56xl0"—  and 
mm  m 


E  =  yjEx^  +  Ey^  +  =  78.7  X 10^  —  (exactly  as  seen  in  Figure  85). 

m 


In  order  to  define  the  effective  distance  between  the  plates; 


,  applied  voltage  15F 
d  = — - ^ — —  = - —  =  lylpm 


electric  field  78  7x10' 


V 


m 


It  is  important  to  note  that  the  effective  distance  is  smaller  in  magnitude  from  the 
real  distance  between  the  plates  where  the  voltage  is  applied. 


Then,  using  Sample  #9,  in  the  configuration  of  Figure  77,  the  electron 
beam  is  applied  to  a  specific  point  on  the  surface  of  the  semiconductor  material. 
A  series  of  images  of  the  intensity  are  taken  (one  with  zero  and  one  with  applied 
voltage),  incrementing  the  applied  voltage  each  time  in  0.5  V  step.  By  evaluating 
the  intensity  ratio  of  the  peak  values  as  described  in  previous  sections  to 
determine  the  local  magnitude  of  the  electric  field  at  the  point  of  charge 
generation,  the  resultant  values  are  obtained  (Figures  86  and  87). 


Normalized  Intensity  Ratio  vs  Applied  Voltage  (Experimental  Data) 
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Figure  86.  Intensity  ratio  vs.  applied  voltage  for  Case  1. 
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From  Figure  87  it  can  be  seen  clearly  that,  as  expected,  the  actual- 
experimental  values  of  the  electric  field  (circles)  in  the  semiconductor  material 
are  smaller  than  the  theoretical  values  (solid  line)  which  are  derived  from  the 
division  of  the  applied  voltage  over  the  effective  distance  between  the  plates.  The 
difference  can  be  taken  as  a  measure  of  the  voltage  drop  at  the  contacts.  From 
Figure  87  it  can  be  seen  that  the  “turn  on”  voltage  of  the  forward  biased  Schottky 
contact  is  at  «  2-3  V  DC  bias.  Below  that  voltage,  very  little  applied  field  exists  in 
the  center  region  of  the  device. 

2.  Case  2 

For  this  case  the  configuration  is  shown  in  Figure  88.  15V  DC  has  been 
applied  between  the  inner  plates  whose  distance  is  185  pm.  The  area  of  interest 
is  illustrated  in  Figure  89  where  a  plot  of  electric  potential  has  been  created  by 
the  use  of  COMSOL  software.  The  color  interpretation  is  given  in  the  color  scale 
shown  at  the  right  of  Figure  89.  The  contour  lines  represent  the  equipotential 
areas.  The  red  circle  shows  the  position  of  the  electron  beam  (not  to  scale). 
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Figure  88.  Configuration  of  Sample  #9  for  defining  the  drop  of  voltage 
between  the  contacts  with  2-point  measurement. 


Figure  89.  COMSOL  simulation  of  the  electric  field  potential  for  Case  2. 
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The  X,  y  and  z  components  of  the  generated  electric  field  between  the 
metal  contacts  which  are  superjacent  to  the  semiconductor  material,  can  be  seen 
in  Figures  90,  91  and  92. 


Figure  90.  COMSOL  simulation  of  the  electric  field  (x-component)  for 

Case  2. 


15V 

Figure  91.  COMSOL  simulation 


•  OV 


the  electric  field  (y-component)  for 
Case  2. 
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Figure  92.  COMSOL  simulation  of  the  electric  field  (z-component)  for 

Case  2. 


By  extracting  the  data  of  Figures  90  through  92  in  the  MATLAB  code  of 
Appendix  C,  graphs  of  the  electric  field  across  a  line  perpendicular  to  the  inner 
plates  that  passes  through  the  position  of  the  electron  beam  (red  circle)  can  be 
seen  in  Figures  93,  94  and  95. 
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Figure  93.  Electric  field  (x-component)  for  Case  2. 
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The  total  electric  field  can  be  seen  in  Figure  96. 
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Figure  96.  Total  electric  fieid  for  Case  2. 


By  determining  the  value  of  the  electric  field  in  the  middle  of  the  distance 
between  the  plates  (i.e.,  185/2=92.5  pm),  in  each  direction,  the  resultant  value  is; 

=  7.827x10"—  ,  Ey  =  lA—,  .Ez  =  7.821x10"  —  and 
mm  m 

E  =  ^JEx^  +  Ey^  =111x10^—  (exactly  as  seen  in  Figure  85). 

m 

In  order  to  define  the  effective  distance  between  the  plates: 

,  applied  voltage  151^ 
d  = - = - =—  =  136//m 

electric  field  1 1 1  x  1 0^  — 

m 

Again  it  should  be  noted  that  the  effective  distance  is  smaller  in  magnitude  from 
the  real  distance  between  the  plates  where  the  voltage  is  applied. 

Next,  by  using  the  specimen  (Sample  #9),  in  the  configuration  of  Figure 
88,  the  electron  beam  is  applied  to  a  specific  point  on  the  surface  of  the 
semiconductor  material.  A  series  of  images  of  the  intensity  are  taken  (one  with 
zero  and  one  with  applied  voltage),  incrementing  the  applied  voltage  each  time  in 
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0.5  V  step.  By  evaluating  the  intensity  ratio  of  the  peak  values  as  described  in 
previous  sections  to  determine  the  local  magnitude  of  the  electric  field  at  the 
point  of  charge  generation,  the  resultant  values  are  obtained  (Figures  97,  98  and 
99). 
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Normalized  Intensity  Ratio  vs  Applied  Voltage  (Experimental  Data) 
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Figure  97.  Intensity  ratio  vs.  applied  voltage  for  Case  2. 
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Figure  99.  Electric  field  vs.  applied  voltage  for  Case  2  (log-log  scale). 

Similarly  to  Case  1  of  the  current  chapter,  the  experimental  values  are 
lower  in  magnitude  than  the  theoretical  values  and  the  “turn  on”  voltage  for  the 
forward  biased  Schottky  diode  is  at  about  2  V  DC  bias. 

B.  COMPARISON  OF  ELECTRIC  FIELD  MAGNITUDE  DETERMINED 

EXPERIMENTALLY  TO  COMSOL  MODELING 

In  this  part  of  the  thesis  it  is  important  to  evaluate  the  accuracy  of  the 
experimental  approach  of  determining  locally  the  magnitude  of  the  electric  field  in 
a  planar  semiconductor  structure,  as  explained  in  Chapter  IV.  One  way  to  do  so 
is  by  comparing  the  results  with  a  simulation  of  an  idealized  structure.  The  latter 
can  be  efficiently  simulated  by  the  use  of  COMSOL  software. 

The  case  study  is  illustrated  in  Figure  100  below  (this  case  has  been 
analyzed  in  Chapter  IV,  paragraph  A-2).  15  V  DC  has  been  applied  between  the 
contacts  and  36  images  taken  at  positions  determined  with  the  red  circles.  In 
addition,  by  using  COMSOL,  the  electric  field  across  a  line  starting  from  position 
1  (PI)  and  ending  at  position  36  (P36)  has  been  evaluated.  The  resultant 
comparison  is  shown  in  Figure  101. 
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Figure  100.  Part  of  Sample  #9  for  comparison  of  the  COMSOL  simulation 

with  the  experimental  method. 


Figure  101.  Comparison  of  the  electric  field  as  a  function  of  position  for  15 

V  DC  applied  bias. 
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From  Figure  101  it  can  be  seen  that  the  experimental  method  used  in  the 
current  work  is  in  good  conditions  with  the  simulations  (about  10%)  and  it  gives 
reasonable  values  of  the  magnitude  of  the  electric  field  locally  based  on  the  real 
characteristics  of  the  material  used.  Local  defects  can  also  play  a  significant  role 
in  decreasing  the  field  at  any  given  spot  and  this  is  something  that  cannot  be 
simulated  in  COMSOL.  Also  it  should  be  noted  that  the  simulation  in  COMSOL 
was  based  on  a  structure  that  has  contacts  made  of  Copper  superjacent  to 
Silicon  semiconductor  material.  The  fluctuations  that  can  be  observed  may  are 
also  due  to  a  small  difference  in  the  position  of  the  line  used  for  COMSOL 
calculations  and  the  position  of  the  electron  beam  during  the  experimental 
approach. 

C.  EVALUATION  OF  THE  ELECTRIC  FIELD  MAGNITUDE  IN 
SEMICONDUCTOR  PLANAR  STRUCTURES  OF  CONSTRICTED 
GEOMETRY 

In  order  to  better  demonstrate  the  potential  capabilities  of  the  experimental 
method  for  determining  the  electric  field  magnitude  as  described  in  this  thesis, 
the  case  illustrated  in  Figure  102  is  explored.  The  motion  of  charge  from  a  region 
of  1  dimensional  (ID)  confinement  to  a  2D  region  is  of  high  interest  in 
preparation  for  the  future  study  of  contacted  nanowires. 


Figure  102.  Part  of  Sample  #9  with  special  geometry  characteristics. 
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In  Figure  102  the  two  upper-right  plates  have  been  bridged  together  and 
the  potential  difference  is  applied  at  points  illustrated  with  the  yellow  lines.  The 
electron  beam  is  applied  at  the  point  specified  with  the  red  circle  and  three 
images  (Figure  103)  are  acquired  with  different  applied  voltages  (zero  bias,  +7  V 
DC  and  -7  V  DC).  The  electron  beam  energy  was  30  keV,  the  probe  current  was 
6E-10  A  and  the  SEM  was  operated  in  spot  mode. 


Figure  103.  Spot  mode  images  at  point  of  charge  generation  with  zero  bias 

(a),  +7  V  DC  bias  (b)  and  -7  V  DC  bias  (c). 


By  using  the  SEM  in  line  mode  (Figure  104)  a  series  of  images  have  been 
acquired  as  demonstrated  in  Figures  105  and  106.  Figure  105  corresponds  to  the 
vertical  red  line  in  Figure  104  and  Figure  106  to  the  green  horizontal  line.  The 
reason  for  using  the  line  mode  is  that  it  gives  more  general  spatial  information  for 
the  behavior  of  the  electric  field  in  the  area  of  interest. 


Figure  104.  SEM  operation  in  line  mode  for  Sample  #9  for  the  case  study. 
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Figure  105.  Line  mode  images  (vertical  red  line)  with  zero  bias  (a),  +7  V  DC 

bias  (b)  and  -7  V  DC  bias  (c). 


(a^  (0  bias) 

(b)  (+7VDC) 

Ip 

(0)  (-7  V  DC) 

(250  pm  X  202  Min) 

(250  Mm  X  202  Mm) 

(250  Mm  X  202  Mm) 

Figure  106.  Line  mode  images  (horizontal  green  line)  with  zero  bias  (a),  +7 

V  DC  bias  (b)  and  -7  V  DC  bias  (c). 


By  the  above  qualitative  analysis  it  has  been  demonstrated  how  important 
a  role  the  geometry  plays  on  the  electric  field  behavior.  Based  on  the  polarity  of 
the  applied  voltage,  the  field  is  either  forced  to  “tunnel”  through  the  space 
between  the  contacts  or  diverges  in  the  opposite  direction.  It  is  important  to  note 
that  the  distance  between  the  “bridged”  contacts  is  about  4  pm  which  is  almost 
equal  to  the  diffusion  length  of  the  minority  carriers. 

In  order  to  quantitatively  determine  the  electric  field  magnitude  at  specific 
points,  the  experimental  method  described  in  previous  sections  is  used.  The 
results  can  be  seen  in  Figures  107  and  108.  In  addition,  in  Figure  108  a 
comparison  between  the  experimental  data  and  a  corresponding  simulation  in 
COMSOL  has  been  done. 
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Figure  107.  E-Field  vs.  distance  for  points  along  (in  the  direction  shown)  of 
the  red  vertical  line  of  Figure  104  with  -  7  V  DC  applied  bias. 


Figure  108.  E-Field  vs.  distance  for  points  along  (in  the  direction  shown)  of 
the  green  horizontal  line  of  Figure  104  with  -  7  V  DC  applied 

bias. 
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By  combining  the  data  of  Figure  107  with  the  qualitative  approach 
demonstrated  in  Figure  105-(c),  it  can  be  seen  that  the  electric  field  magnitude  at 
points  constrained  between  the  two  “bridged”  contacts  is  much  lower  than  the 
values  of  the  field  in  the  “open”  area  between  the  points  where  the  bias  is 
applied.  In  Figure  107  the  first  10  field  values  correspond  to  points  inside  the 
constrained  area  between  the  “bridged”  contacts  and  the  remaining  15  values 
correspond  to  points  in  the  2D  area.  The  experimental  approach  shows  existence 
of  E-field  in  places  where  COMSOL  MULTIPHYSICS  did  not. 

For  the  data  of  Figure  108,  combined  with  that  of  Figure  106-(c)  it  can  be 
seen  that  the  lowest  value  of  the  electric  field  is  exactly  at  the  point  where  the 
constrained  area  between  the  “bridged”  contacts  is  connected  with  the  “open” 
area.  All  the  above  areas  are  illustrated  in  Figure  109.  The  experimental  data 
have  also  been  compared  with  the  results  of  a  simulation  made  by  COMSOL 
MULTIPHYSICS  in  a  similar  semiconductor  structure.  The  trend  of  the 
experimental  data  is  consistent  with  the  simulation  but  the  magnitude  has  been 
scaled  in  order  to  include  the  effects  of  fluctuations  in  contact  resistance,  an 
effect  that  cannot  be  simulated. 


Figure  109.  Areas  of  interest  for  determination  of  the  electric  field 

behavior. 

This  last  experiment  reveals  useful  details  associated  with  charge  injection 
in  confined  regions,  details  that  cannot  be  otherwise  measured  electrically. 
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VI.  CONCLUSIONS  AND  SUGGESTIONS  FOR  FURTHER 

RESEARCH 

A.  SUMMARY  AND  CONCLUSION 

Photonic  devices  can  be  divided  in  the  following  categories  with  respect  to 
their  principles  of  operation:  Light-emitting  diodes  (LEDs)  and  lasers  (light 
amplification  by  stimulated  emission  of  radiation)  convert  electrical  energy  to 
optical  energy.  Photodetectors  detect  optical  signals  and  solar  cells  convert 
optical  energy  into  electrical  energy.  Specifically  the  operation  of  photodetectors, 
which  are  used  extensively  in  various  kinds  of  optical  sensors,  is  based  on  carrier 
generation  by  incident  light,  carrier  transport  and  interaction  of  current  with  the 
external  circuit  to  provide  the  output  signal  [12].  A  wide  range  of  devices  depend 
upon  the  fabrication  of  contacts  that  can  uniformly  inject  and  remove  charge.  The 
understanding  of  local  transport  behavior  is  essential  for  the  development  of  new 
devices  in  all  these  areas. 

The  work  described  in  this  thesis  is  based  on  cathodoluminescence 
scanning  electron  microscopy.  A  technique  for  imaging  the  transport  of  free 
charge  in  planar  semiconductor  structures  has  been  developed  in  our  laboratory. 
This  technique  has  been  used  to  image  the  motion  of  minority  carriers  in  a 
heavily  doped  p-type  AIGaAs/GaAs/AIGaAs  heterostructure,  under  a  range  of 
applied  electric  fields.  By  generating  charge  at  a  fixed  point,  it  is  possible  to 
monitor  the  subsequent  drift  and  diffusion  by  imaging  the  distribution  of  the  free 
carrier  recombination  path  as  illustrated  in  Figure  110  below.  In  this  figure,  the 
SEM  is  operated  in  spot  mode  and  three  images  have  been  acquired,  with  zero, 
DC  and  AC  bias  in  Figures  (a),  (b)  and  (c)  respectively.  In  the  following  3D 
graphs,  the  intensity  of  light  at  the  point  of  generation  is  acquired  for  each  case 
mentioned  above.  Additionally,  the  direction  of  the  drift  length  due  to  the  applied 
bias  is  efficiently  illustrated. 
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Figure  110.  SEM  images  in  spot  mode  operation  for  point  charge 

generation. 


The  technique  described  in  this  thesis  allows  the  mapping  of  local  electric 
fields  with  high  resolution  by  analyzing  the  transport  behavior  under  the 
combined  influence  of  both  diffusion  and  drift  predicted  by  modeling  the  drift  and 
diffusion  in  2D  following  generation  at  a  point  source.  Subsequent  measurements 
of  electric  field  profiles  have  been  performed  in  both  uniform  and  non-uniform 
regions.  The  acquired  results  agree  (to  within  about  10%)  with  simulations  made 
in  COMSOL  MULTIPHYSICS  software  for  similar  geometry  semiconductor 
planar  structures.  Fluctuations  between  the  experimental  and  simulated  method 
are  due  to  material  variations,  defects  and  actual  characteristics  of  the  material 
used  (e.g.,  type  and  quantity  of  doping),  conditions  that  cannot  be  easily 
simulated. 

Additionally,  the  experimental  method  described  in  the  thesis  allows  the 
accurate  measurement  of  the  contact  resistance  with  only  a  two-point  contact. 
The  technique  illustrates  the  device’s  characteristics  by  determining  the  exact 
activation  point  of  the  forward  biased  Schottky  diode  and  further  deviations  from 
an  ideal  linear  l-V  behavior. 
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Finally,  the  method  was  demonstrated  for  a  geometry  where  the  non¬ 
uniformity  of  the  electric  field  emphasizes  the  necessity  of  having  an 
experimental  method  in  order  to  determine  fields  with  high  resolution.  The 
technique  described  in  this  thesis  is  a  powerful  “tool”  for  testing  the  reliability  of 
electronic  devices  in  a  radiative  environment  or  in  conditions  that  may  alter  their 
electrical  behavior. 

B.  SUGGESTIONS  FOR  FURTHER  RESEARCH 

1 .  Reliability  of  Semiconductor  Devices 

After  the  industrial  processing  of  a  wafer  containing  a  large  number  of 
integrated  circuits,  the  term  “yield”  refers  to  the  fraction  of  successful  products 
that  survive  final  screening  and  meet  specifications  relative  to  the  total  number 
that  started  out  to  be  processed  or  manufactured.  The  fact  that  yielded  products 
have  passed  the  last  screen  does  not  mean  they  do  not  contain  potential  or 
latent  defects  [13].  The  reliability  of  these  devices  is  directly  based  on  the 
existence  of  defects  that  influence  their  electrical  behavior. 

The  experimental  technique  that  is  described  in  this  thesis  allows  for 
detecting  these  defects  not  only  optically  (by  using  the  corresponding  CCD 
camera)  but  also  quantitatively  by  determining  the  magnitude  of  the  electric  field 
with  high  resolution.  Since  defects  influence  the  electrical  behavior  of 
semiconductor  devices,  the  capability  of  mapping  the  electric  field  locally  will 
detect  the  potential  existence  of  defects. 

Additional  to  the  existence  of  defects,  a  phenomenon  that  has  a  significant 
effect  on  the  reliability  of  semiconductor  devices  is  electromigration. 
Electromigration  is  characterized  by  the  migration  of  metal  atoms  in  a  conductor 
through  which  large  direct-current  densities  pass.  Passage  of  high  current 
densities  through  interconnects  causes  time-dependent  mass-transport  effects 
that  are  manifested  as  surface  morphological  changes. 
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2.  Electromigration  Issues 

At  a  fundamental  level,  electromigration  involves  the  interaction  between 
current  carriers  and  migrating  atoms.  Since  any  metal  that  leaves  the 
semiconductor  interface  cannot  be  replenished  [13],  electromigration  at  contacts 
is  of  great  importance.  During  the  experimental  work  done  for  this  thesis,  it  has 
been  noticed  that  after  a  period,  there  was  a  significant  drop  in  the  resistance  of 
the  specimen  used  (Sample  #9)  resulting  in  the  degradation  of  its  performance. 
In  order  to  understand  the  reasons  for  this  material  degradation,  a  series  of 
experiments  has  been  done  in  order  to  force  electromigration  to  happen  in 
Sample  #9.  A  characteristic  example  of  the  experimental  setup  for  this  purpose  is 
illustrated  in  Figure  111. 


Figure  111.  Sample  #9  for  electromigration  study  (each  image  is  874  pm  X 

589  pm). 


The  electron  beam  is  positioned  as  shown  in  Figure  111  (a)  in  the  non- 
uniform  region  of  the  electric  field  created  by  the  application  of  12.5  V  DC  bias 
between  the  contacts.  Electron  beam  energy  was  30  keV  and  the  probe  current 
was  3X10'^°  A.  The  specimen  remained  under  the  radiation  of  the  electron  beam 
for  2h  30min  without  changing  the  position  of  either  the  beam  or  the  sample.  A 
series  of  images  have  been  taken  in  periodic  intervals  in  order  to  detect  possible 
deformations  of  the  contacts.  Part  (b)  of  Figure  111  represents  the  final  image  of 
the  specimen  after  2h  30min  of  radiation  (taken  without  applying  bias).  By 
analyzing  the  images  in  high  resolution,  the  only  difference  between  parts  (a)  and 
(b)  is  the  surface  field  distribution  which  disappears  when  there  is  no  bias. 
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Naturally,  after  2h  30min  of  radiation  there  was  a  significant  increase  in  the 
temperature  of  the  specimen  (from  room  temperature  to  almost  77°C).  It  is  very 
interesting  at  this  point  to  notice  the  difference  of  the  electric  field  magnitude  due 
to  temperature  increase  with  respect  to  time  (Figure  112). 


Figure  112.  E-Field  values  with  respect  to  time  due  to  temperature 

increase. 

In  Figure  1 12  it  can  be  seen  that  the  electric  field  within  the  semiconductor 
material  (at  the  point  of  the  incident  electron  beam)  increases  as  time  passes 
due  to  the  increase  in  temperature.  That  means  that  the  contact  resistance 
decreases  due  to  heating.  Another  point  of  interest  is  the  existence  of  some 
white  spots  in  the  perimeter  of  the  contact  where  the  bias  was  applied  [Figure 
111  (a)].  These  white  spots  are  points  where  injected  charge  density  is  extremely 
high  and  represent  potential  points  where  electromigration  can  be  accelerated. 
This  can  be  seen  clearly  in  Figure  113  below  where  with  the  application  of  bias, 
high  charge  density  is  noticed  in  the  corner  of  the  plate  (demonstrated  with  the 
yellow  arrow). 
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Figure  113.  Points  with  high  injected  charge  density  (each  image  is  247 

|jm  X  215  |jm). 

The  reasons  that  led  to  the  degradation  of  Sample  #9  may  not  be  related 
only  with  the  electromigration  mechanisms.  Possible  degradation  of  the  Schottky 
contacts  used  in  the  specimen  (Au/Ti/Alo.sGao.sAs/GaAs  capping  layer)  may 
change  the  electrical  behavior  of  the  material  due  to  changes  to  the  metal- 
semiconductor  interface.  Additionally,  the  three  kind  of  tests  described  in  [13] 
concerning  the  electromigration: 

•  the  temperature-ramp  resistance  analysis  to  characterize 

electromigration  (TRACE), 

•  the  breakdown  energy  of  metal  (BEM)  and 

•  the  standard  wafer  -  level  electromigration  acceleration  test 
(SWEAT),  all  are  techniques  referred  especially  to  the  wafer  level, 
requiring  simultaneously  specific  laboratory  conditions  and 
capabilities. 

Since  in  electromigration,  metal  atoms  are  forced  to  migrate  by  a  large 
current  density  and  such  areas  with  high  charge  density  can  be  identified  by  the 
SEM  imaging  capabilities  as  shown  in  Figures  111  and  113  above,  an  effective 
method  of  studying  the  electromigration  issues  can  be  developed.  It  is  also 
essential  to  study  the  effects  of  the  electromigration  first  in  small  cross-sectional 
interconnects  where  high  current  densities  are  expected.  The  collision  between 
fast  moving  electrons  having  high  drift  velocity  with  the  metal  ions  can  result  in 
material  depletion  (voids),  material  accumulation  (hillocks)  or  in  an  interconnect 
failure  [3]. 
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Finally,  another  case  has  been  studied  in  order  to  explore  the  effect  of 
“cycling”  the  applied  bias.  The  term  “cycle”  means  that  bias  is  applied,  then 
switched  off  and  then  re-applied.  The  structure  used  is  illustrated  in  Figure  114 
below.  Thermal  excitations  of  the  lattice  structure  due  to  abrupt  voltage 
differences  are  anticipated  to  influence  the  electrical  behavior  of  the  material. 


Figure  114.  Part  of  Sample  #9  for  studying  the  effects  of  thermal  excitation 

(874  pm  X  589  pm) 

In  Figure  114  above,  the  position  of  the  electron  beam  is  shown  with  the 
red  circle.  Electron  beam  energy  is  at  30  keV  and  the  probe  current  at  3E-9  A. 
The  applied  voltage  was  -15  V  DC  bias  and  it  was  switched  on  and  off  for  30 
times  continuously  with  15  sec  intervals  between  each  cycle.  The  acquired 
results  are  shown  in  Figure  115.  It  can  be  seen  from  the  figure  that  the  values  of 
the  electric  field  is  not  in  an  horizontal  line,  as  expected  in  a  theoretical  situation 
like  this,  but  they  have  some  fluctuations  («  4%)  due  to  thermal  excitation.  By 
microscopic  analysis,  there  were  not  any  signs  of  electromigration  or  changes  in 
the  formation  of  the  structure. 
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Figure  115.  E-Field  values  vs.  times  of  power  cycling. 


In  summary,  this  work  has  successfully  demonstrated  that  transport 
imaging  can  be  used  as  a  quantitative  method  for  electric  field  mapping.  The 
technique  has  potential  in  the  study  of  defect  behavior  and  electromigration, 
areas  with  extreme  scientific  significance  in  the  micro-electronics  industry. 


104 


APPENDIX  A.  MATLAB  CODE  FOR  EVALUATING  THE  DC 

MATHEMATICAL  MODEL 


%FIND  THEORETICAL  DC  MODEL  EQUATION 
%CREATED  BY  PAVLOS  ANDRIKOPOULOS 

%READ  EXCEL  FILES  AND  EXTRACT  DATA 
A=xlsread( 'MODEL_DATA' ) 
xl=  A(:,l) 

[mxl,  r] =max (xl) 
nunxl=max  (mxl ) 
x2=  A ( : , 2 ) 

[mx2,  r] =max (x2) 
nunx2=max  (mx2 ) 
x3=  A ( : , 3 ) 

[mx3,  r] =max (x3) 
nunx3=max  (mx3) 
x4=  A ( : , 4) 

[mx4,  r]  =max  (x4) 
nunx4=max  (mx4 ) 
x5=  A ( : , 5 ) 

[mx5,  r] =max (x5) 
nunx5=max  (mx5) 
x6=  A ( : , 6 ) 

[mx6, r] =max (x6) 
nunx6=max  (mx6) 
x7=  A(:,7) 

[mx7,  r]  =max  (x7) 
nunx7=max  (mx7 ) 
x8=  A ( : , 8 ) 

[mx8,  r]  =max  (x8) 
nunx8=max  (mx8 ) 
x9=  A(  :  ,  9) 

[mx9, r] =max (x9) 
nunx9=max  (mx9) 
xlO=  A( : , 10) 

[mxlO, r] =max (xlO) 
nunxlO=max  (mxlO ) 
xll=  A ( : , 11) 

[mxll, r] =max (xll) 
mmxl l=max (mxll ) 
xl2=  A ( : , 12) 

[mxl2, r] =max (xl2) 
mmxl2=max (mxl2 ) 
xl3=  A( : , 13) 

[mxl3, r] =max (xl3) 
mmxl3=max (mxl3) 
xl4=  A ( : , 14) 

[mxl4, r] =max (xl4) 
mmxl4=max (mxl 4 ) 
xl5=  A( : , 15) 

[mxl5, r] =max (xl5) 
mmxl5=max (mxl5) 
xl6=  A( : , 16) 
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[mxl6, r] =max (xl6) 
nunxl  6=max  (mxl  6 ) 
xl7=  A ( : , 17) 
[mxl7,  r] =max (xl7) 
nunxl7=max  (mxl  7 ) 
xl8=  A( : , 18) 
[mxl8,  r] =max (xl8) 
nunxl8=max  (mxl  8 ) 
xl9=  A( : , 19) 
[mxl9, r] =max (xl9) 
mmxl 9=max (mxl 9 ) 
x20=  A ( : , 20 ) 
[mx20,  r] =max (x20) 
mmx20=max (mx20 ) 
x21=  A ( : , 21) 
[mx21, r] =max (x21) 
mmx2 l=max (mx2 1 ) 
x22=  A ( : , 22) 
[mx22, r] =max (x22) 
mmx22=max (mx22 ) 
x23=  A ( : , 23) 
[mx23, r] =max (x23) 
mmx2  3=max (mx2  3 ) 
x24=  A ( : , 24) 
[mx24,  r] =max (x24) 
mmx24=max (mx24 ) 
x25=  A ( : , 25) 
[mx25, r] =max (x25) 
mmx2  5=max (mx2  5 ) 
x2  6=  A ( : , 2  6 ) 
[mx26, r] =max (x26) 
mmx2  6=max (mx2  6 ) 
x27=  A ( : , 27) 
[mx27,  r] =max (x27) 
mmx27=max (mx27 ) 
x28=  A ( : , 28 ) 
[mx28,  r] =max (x28) 
mmx28=max (mx28 ) 
x29=  A(:,29) 
[mx29, r] =max (x29) 
mmx2  9=max (mx2  9 ) 
x30=  A ( : , 30 ) 
[mx30,  r] =max (x30) 
mmx30=max (mx30 ) 
x31=  A( : , 31) 
[mx31, r] =max (x31) 
mmx31=max (mx31 ) 
x32=  A ( : , 32 ) 
[mx32, r] =max (x32) 
mmx32=max (mx32 ) 
x33=  A ( : , 33 ) 
[mx33, r] =max (x33) 
mmx33=max (mx33) 
x34=  A( : , 34) 
[mx34,  r] =max (x34) 
mmx34=max (mx34 ) 
x35=  A ( : , 35 ) 
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[mx35, r] =max (x35) 
nunx35=max  (mx35) 
x3  6=  A ( : , 3  6 ) 

[mx36, r] =max (x36) 
nunx36=max  (mx36) 

%CREATE  MATRIX  FOR  MAX  VALUES 

X=[mmxl  mmx2  mmx3  mmx4  mmx5  mmx6  mmx7  mmx8  mmx9  mmxlO  mmxll  mmxl2  mmxl3 
nunxl4  mmxl5  mmxl6  mmxl7  mmxl8  mmxl9  mmx20  mmx21  mmx22  mmx23  mmx24  mmx25 
nunx26  mmx27  mmx28  mmx29  mmx30  mmx31  mmx32  mmx33  mmx34  mmx35  mmx36] 

%DETERMINE  THE  X  VECTOR 
x=(-30:0.5:30) 

%CREATE  THE  GRAPH  INTENSITY  VS  x 
plot(x,A(:,l), 'b') 
hold  on 

plot (x,A( : , 4) , 'g' ) 
hold  on 

plot (x, A ( : , 6) , ' r ' ) 
hold  on 

plot (x, A ( : , 9) ,  ' c  '  ) 
hold  on 

plot(x,A(:,ll), 'm') 
hold  on 

plot (x,A( : , 13) , 'y' ) 
hold  on 

plot (x,A( : , 15) , 'k' ) 
hold  on 

plot (x,A( : , 17) ,  'b-'  ) 
hold  on 

plot (x,A( : , 19) , 'g-' ) 
hold  on 

plot (x, A ( : , 22) , ' r- ' ) 
hold  on 

plot(x,A(:,25), 'c-') 
hold  on 

plot (x, A ( : , 27 ) , 'm- ' ) 
hold  on 

plot (x, A ( : , 29) , ' y- ' ) 
hold  on 

plot (x,A( : , 32) , 'k-' ) 
hold  on 

plot (x,A( : , 34) ,  'b.  '  ) 
hold  on 

plot (x,A( : , 36) , 'g. ' ) 
hold  off 

grid  on 
xlabel ( ' X ' ) 

ylabel (' Intensity  Ratio') 

title (' Normalized  Intensity  Ratio  (Theoretical  Model) ') 

legend ( ' 0  V/\mum' , ' 30E-4  V/\mum' , ' 44E-4  V/\mum' , ' 70E-4  V/\mum' , ' lOOE-4 
V/\mum', '131.50E-4  V/\mum', '170E-4  V/\mum', '209.50E-4  V/\mum', '250E-4 
V/\mum',  '313.50E-4  V/\mum',  '370E-4  V/\mum',  '412E-4  V/\mum',  '450E-4 
V/ \mum ' , ' 52 lE-4  V/\mum' , ' 560E-4  V/\mum' , ' 600E-4  V/\mum',2) 
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%CURVE  FITTING  OF  THF  RESULTS 
B=xlsread( ' VOLTAGE_DATA ' ) 

Eold=  B  (  : , I) 

E=Eold' 

Iold=X 

Inorm=Iold/0 .0I7393I86II4I5I40 
n=2 ; 

p=polyf it (Inorm, E, n) 
xi= (0 . 40:0.01:1.0); 
yi=polyvaI (p, xi) ; 
figure 

plot ( Inorm, E, 'o',xi,yi, 
grid  on 

xlabel ( ' Peak  Intensity  Ratio '), ylabel (' E-Field  (V/\mum)') 

title ('Curve  Fit  for  DC  Theoretical  Model') 

text  (0 . 45,  0 . 005,  '  (0.0842*I''2)-(0.2162*I)  +  (0.1351)  ') 

legend (' Theoretical  Model ',' Curve  Fitting', 1) 

figure 

plot  (E, Inorm, 'o') 
grid  on 

xlabel (' E-Field  (V/ \mum) '), ylabel (' Peak  Intensity  Ratio') 
title  (' Intensity  ratio  vs  Applied  Electric  Field') 

%  Etheory= (0 . 0842*1. *I)-(0. 21 62*1) +0.1351 
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APPENDIX  B.  MATLAB  CODE  FOR  MAPPING  ELECTRIC  FIELD 
VIA  TRANSPORT  IMAGING  TECHNIQUE 


%MATLAB  CODE  FOR  FINDING  THE  ELECTRIC  FIELD  FROM  SEM  IMAGES 
%DC  CURRENT 

%CREATED  BY  ANDRIKOPOULOS  PAVLOS 

%POSITION_I-OV  UP  TO  POSITION_23-OV 
%READ  TIE  FILES  AND  EXTRACT  THE  MAX  VALUE 

ZI=imread ( 'P1_0V.TIF ' ) 

ZI=double (ZI) 

[YI,X] =max (ZI) 
nunxl=max  ( YI ) 

Z2=imread ( 'P2_0V.TIF ' ) 

Z2=double ( Z2 ) 

[ Y2 , X] =max ( Z2 ) 
nunx2=max  ( Y2 ) 

Z3=imread ( 'P3_0V.TIF ' ) 

Z3=double ( Z3 ) 

[Y3,  X] =max (Z3) 
nunx3=max  ( Y3 ) 

Z4=imread ( 'P4_0V.TIF ' ) 

Z4=double ( Z4 ) 

[Y4,  X] =max (Z4) 
nunx4=max  ( Y4 ) 

Z5=imread ( 'P5_0V.TIF ' ) 

Z5=double ( Z5 ) 

[Y5,  X] =max (Z5) 
nunx5=max  ( Y5 ) 

Z6=imread ( 'P6_0V.TIF ' ) 

Z6=double (Z6) 

[Y6, X] =max (Z6) 
nunx6=max  (Y6) 

Z7=imread ( 'P7_0V.TIF ' ) 

Z7=double ( Z7 ) 

[ Y7 , X] =max  (ZI) 
nunx7=max  ( Y7 ) 

Z8=imread ( 'P8_0V.TIF ' ) 

Z8=double ( Z8 ) 

[ Y8 , X] =max ( Z8 ) 
nunx8=max  ( Y8 ) 

Z9=imread ( 'P9_0V.TIF ' ) 

Z9=double (Z9) 

[Y9, X] =max (Z9) 
nunx9=max  ( Y9 ) 

ZIO=imread( 'PI0_0V.TIF' ) 

ZI0=double  (ZIO) 

[YI0,X] =max (ZIO) 
nunxIO=max  (YIO ) 

ZII=imread( 'P11_0V.TIF' ) 

ZII=double (ZII) 

[YII,X] =max (ZII) 
nunxi  I=max  ( YI  I ) 


109 


Z12=imread( 'P12_0V.TIF' ) 

Z12=double (Z12) 

[Y12,X] =max (Z12) 
nunxl2=max  ( Y12 ) 

Z13=imread( 'P13_0V.TIF' ) 

Z13=double (Z13) 

[Y13,X] =max (Z13) 
nunxl3=max  (Y13) 

Z14=imread( 'P14_0V.TIF' ) 

Z14=double  (Z14) 

[Y14,X] =max (Z14) 
nunxl4=max  (Y14 ) 

Z15=imread( 'P15_0V.TIF' ) 

Z15=double (Z15) 

[Y15,  X] =max (Z15) 
nunxl5=max  (Y15) 

Z16=imread( 'P16_0V.TIF' ) 

Z16=double  (Z16) 

[Y16, X] =max (Z16) 
nunxl6=max  (Y16) 

Z17=imread( 'P17_0V.TIF' ) 

Z17=double  (Z17) 

[ Y17 , X] =max ( Z17 ) 
nunxl7=max  (Y17 ) 

Z18=imread( 'P18_0V.TIF' ) 

Z18=double  (Z18) 

[Y18,X] =max (Z18) 
nunxl8=max  (Y18 ) 

Z19=imread( 'P19_0V.TIF' ) 

Z19=double  (Z19) 

[Y19,X] =max (Z19) 
nunxl9=max  (Y19) 

Z20=imread( 'P20_0V.TIF' ) 

Z20=double  (Z20) 

[ Y20 , X] =max ( Z20 ) 
nunx2  0=max  (Y2  0 ) 

Z21=imread( 'P21_0V.TIF' ) 

Z21=double (Z21) 

[Y21,X] =max (Z21) 
nunx2  l=max  ( Y2 1 ) 

Z22=imread( 'P22_0V.TIF' ) 

Z22=double (Z22) 

[  Y22 , X] =max ( Z22 ) 
nunx22=max  ( Y22 ) 

Z23=imread( 'P23_0V.TIF' ) 

Z23=double (Z23) 

[Y23, X] =max (Z23) 
nunx23=max  (Y23) 

%P0SITI0N_1-15V  UP  TO  POSITION_23-15V 
%READ  TIF  FILES  AND  EXTRACT  THE  MAX  VALUE 

Al=imread ( 'P1_15V.TIF ' ) 

Al=double (A1 ) 

[B1 , X] =max (A1 ) 
nunyl=max  (B1 ) 
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A2=imread ( 'P2_15V.TIF '  ) 
A2=double (A2) 

[B2 ,  X] =max  (A2 ) 
nuny2=max  (B2 ) 

A3=imread ( 'P3_15V.TIF ' ) 
A3=double (A3 ) 

[B3, X] =max (A3) 
nuny3=max  (B3 ) 

A4=imread ( 'P4_15V.TIF ' ) 
A4=double (A4) 

[B4 , X] =max (A4 ) 
mmy4=max (B4 ) 

A5=imread ( ' P5_15V. TIF ' ) 
A5=double (A5 ) 

[B5, X] =max (A5) 
nuny5=max  (B5 ) 

A6=imread ( 'P6_15V.TIF ' ) 
A6=double (A6 ) 

[B6, X] =max (A6) 
nuny6=max  (B6 ) 

A7=imread ( 'P7_15V.TIF ' ) 
A7=double (A7) 

[B7 , X] =max (A7 ) 
nuny7=max  (B7 ) 

A8=imread ( 'P8_15V.TIF ' ) 
A8=double (A8 ) 

[B8 , X] =max (A8 ) 
nuny8=max  (B8 ) 

A9=imread ( 'P9_15V.TIF ' ) 
A9=double (A9 ) 

[B9, X] =max (A9) 
nuny9=max  (B9 ) 

A10=imread( 'P10_15V.TIF' ) 
A10=double (AlO) 

[B10,X] =max (AlO) 
nunylO=max  (BIO ) 

All=imread( 'P11_15V.TIF'  ) 
All=double (All) 

[B11,X] =max (All) 
mmyl l=max (B1 1 ) 

A12=imread( 'P12_15V.TIF' ) 
A12=double (A12) 

[B12,X] =max (A12) 
nunyl2=max  (B12 ) 

A13=imread( 'P13_15V.TIF'  ) 
A13=double (A13) 

[B13,X] =max (A13) 
nunyl3=max  (B13) 

A14=imread( 'P14_15V.TIF'  ) 
A14=double (A14) 

[B14,X] =max (A14) 
nunyl4=max  (B14 ) 

A15=imread( 'P15_15V.TIF'  ) 
A15=double (A15) 

[B15, X] =max (A15) 
nunyl5=max  (B15) 

A16=imread( 'P16_15V.TIF' ) 


A16=double (A16) 

[B16, X] =max (A16) 
nunyl6=max  (B16) 

A17=imread( 'P17_15V.TIF'  ) 

A17=double (A17) 

[B17 ,  X] =max (A17 ) 
nunyl7=max  (B17 ) 

A18=imread( 'P18_15V.TIF' ) 

A18=double (A18) 

[B18,X] =max (A18) 
nunyl8=max  (B18 ) 

A19=imread( 'P19_15V.TIF' ) 

A19=double (A19) 

[B19,X] =max (A19) 
nunyl9=max  (B19) 

A20=imread( 'P20_15V.TIF'  ) 

A20=double (A20) 

[B2  0 ,  X] =max (A2  0 ) 
nuny2  0=max  (B2  0 ) 

A21=imread( 'P21_15V.TIF'  ) 

A2 l=double (A2 1 ) 

[B21,X] =max (A21) 
nuny2 l=max  (B2 1 ) 

A22=imread( 'P22_15V.TIF'  ) 

A22=double (A22) 

[B22 ,  X] =max (A22 ) 
nuny22=max  (B22 ) 

A23=imread( 'P23_15V.TIF'  ) 

A23=double (A23) 

[B23, X] =max (A23) 
nuny23=max  (B23) 

%CREATE  MATRIX  FOR  MAX  VALUES  OF  0  VOLTS 

X=[mmxl  mmx2  mmx3  mmx4  mmx5  mmx6  mmx7  mmx8  mmx9  mmxlO  mmxll  mmxl2  mmxl3 
nunxl4  mmxlS  mmxlG  mmxl7  mmxl8  mmxl9  mmx20  mmx21  mmx22  mmx23] 

%CREATE  MATRIX  FOR  MAX  VALUES  OF  15  VOLTS 

Y=[mmyl  mmy2  mmy3  mmy4  mmy5  mmy6  mmy7  mmy8  mmy9  mmylO  mmyll  mmyl2  mmyl3 
nunyl4  mmyl5  mmylG  mmyl7  mmyl8  mmyl9  mmy20  mmy21  mmy22  mmy23] 

%CALCULATE  NORMALIZED  INTENSITY 
Io=240 

1= (Y-Io) . / (X-Io) 

%CALCULATE  VALUES  FOR  DISTANCE 
%SET  VALUES  FOR  X  (H=HORIZONTAL) 

H= [10. 02, 10. 05, 10. 07, 10. 09, 10. 11, 10. 13, 10. 15, 10. 17, 10. 19, 10. 21, 10. 23, 10 
.25,  10. 27,  10. 29,  10. 31,  10. 33,  10. 35,  10. 37,  10. 39,  10. 41,  10. 51,  10. 53,  10. 55] 
%SET  VALUES  FOR  Y  (V=VERTICAL) 

V= [13. 80, 13. 80, 13. 80, 13. 80, 13. 80, 13. 80, 13. 80, 13. 80, 13. 80, 13. 80, 13. 80, 13 
.80,  13. 80,  13. 80,  13. 80,  13. 78,  13. 78,  13. 78,  13. 78,  13. 78,  13. 78,  13. 78,  13. 78] 

%EVALUATE  THE  DISTANCE  FROM  INITIAL  POINT 
D=sqrt ( (H-10.02) .* (H-10 . 02 ) + ( 13 . 80-V) .* (13.80-V) ) 
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%CREATE  THE  GRAPH 
plot (D, I , ' bo ' ) 
grid  on 

xlabel (' Distance  (mm) ') 

ylabel (' Normalized  Intensity  Ratio') 

title (' Normalized  Intensity  Ratio  vs  Distance') 


%EVALUATION  OF  THEORITICAL 
CORRESPONDING  GRAPH 


MODEL  FOR  DC 


AND 


CREATION  OF 


THE 


[0 

0 . 001 

0 . 002 

0 . 003 

0 . 004 

0 . 0044 

0 . 005 

0.006  0. 

007 

008 

0 .01 

0 .012 

0 . 0132 

0 . 015 

0.017 

0 . 02 

0.021  0. 

023 

025 

0 . 027 

0.029 

0 . 0314 

0 . 033 

0 . 035 

0 . 037 

0.039  0.0 

412 

043 

0 . 045 

0 . 047 

0 . 049 

0 . 0521 

0.054  0 

.  056  0 

.058  0.06] 

r 

ld=[0, 

.  0173931861141514 

0 . 017368623726892 

4  0.0 

172957143671 

621 

0171766638980226 

0.017 

0149206611566 

0 . 0 

169392356157 

452 

0 . 0168149071852186 
0 . 0162940163629757 
0 .0150242440677031 
0 . 0131890887893245 
0 .0121406835833828 
0 . 0109547598726452 
0 . 0100454544116619 
0 . 00921432745937045 
0 . 00850034248019468 


0 .0166530028607825 
0 .0158421498452957 
0 .0145185911735649 
0 .0129837833235535 
0 .0117481446599407 
0 .0106754611667708 
0 . 00975404133410867 
0 . 00896439018545242 


0 . 00818920174584799 
0.00783888525452678  0.00767012138997172  0.00750822234977514]; 
Inorm=Iold/0 .01739318611415140 


0 . 0164875036781727 
0 . 0153323559252167 
0 .0139712700358284 
0 . 0125496590763713 
0 .0113730073666232 
0 . 0103523335916252 
0 . 00945037407473868 
0 . 00872663147449506 
0 . 00801488318079925 


n=2 ; 

p=polyf it (Inorm, E, n) 
xi= (0 . 40:0.01:1.0); 
yi=polyval (p, xi) ; 


figure 

plot ( Inorm, E, 'o',xi,yi, '-') 
grid  on 

xlabel (' Intensity  Ratio '), ylabel (' E-Field  (V/mum') 

title ('Curve  Fit  for  DC  Theoretical  Model') 

text (0.42,0.005, 'Etheory= (0 . 0842*I^2)-(0. 21 62*1) +0.1351') 

%APPLY  CORRECTION 

Etheory= (0 . 0842*1. *I)-(0. 21 62*1) +0.1351 


%CREATE  THE  FINAL  GRAPH 
figure 

plot (D, I , 'o') 
grid  on 

xlabel (' Distance  (mm) '), ylabel (' Normalized  Intensity  Ratio') 
title  (  'Experimental  Measurement ' ) 

figure 

plot  (D, Etheory,  'o') 
grid  on 

xlabel (' Distance  (mm) '), ylabel (' E-Field  (V/\mum)') 
title (' E-Field  vs  Distance') 
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APPENDIX  C.  MATLAB  CODE  FOR  EVALUATION  OF  CONTACT 

RESISTANCE 


%EVALUATION  OF  DROP  OF  VOLTAGE  BETWEEN  THE  CONTACT  PLATES 
%DC  CURRENT 

%CREATED  BY  ANDRIKOPOULOS  PAVLOS 

%POSITION  1,  MEASUREMENTS  1-31,  OV  (DC) 

%READ  TIE  FILES  AND  EXTRACT  THE  MAX  VALUE 

Zl=imread ( 'P1_0V.TIF ' ) 

Zl=double (Zl) 

[Y1,X] =max (Zl) 
nunxl=max  ( Y1 ) 

Z2=imread ( 'P2_0V.TIF ' ) 

Z2=double ( Z2 ) 

[ Y2 , X] =max ( Z2 ) 
nunx2=max  ( Y2 ) 

Z3=imread ( 'P3_0V.TIF ' ) 

Z3=double ( Z3 ) 

[Y3,  X] =max (Z3) 
nunx3=max  ( Y3 ) 

Z4=imread ( 'P4_0V.TIF ' ) 

Z4=double ( Z4 ) 

[Y4,  X] =max (Z4) 
nunx4=max  ( Y4 ) 

Z5=imread ( 'P5_0V.TIF ' ) 

Z5=double ( Z5 ) 

[Y5,  X] =max (Z5) 
nunx5=max  ( Y5 ) 

Z6=imread ( 'P6_0V.TIF ' ) 

Z6=double (Z6) 

[Y6, X] =max (Z6) 
nunx6=max  (Y6) 

Z7=imread ( 'P7_0V.TIF ' ) 

Z7=double ( Z7 ) 

[ Y7 , X] =max  (Zl) 
nunx7=max  ( Y7 ) 

Z8=imread ( 'P8_0V.TIF ' ) 

Z8=double ( Z8 ) 

[ Y8 , X] =max ( Z8 ) 
nunx8=max  ( Y8 ) 

Z9=imread ( 'P9_0V.TIF ' ) 

Z9=double (Z9) 

[Y9, X] =max (Z9) 
nunx9=max  ( Y9 ) 

Z10=imread( 'P10_0V.TIF' ) 

Z10=double  (ZIO) 

[Y10,X] =max (ZIO) 
nunxlO=max  (YIO ) 

Zll=imread( 'P11_0V.TIF' ) 

Zll=double (Zll) 

[Y11,X] =max (Zll) 
nunxl  l=max  ( Y1 1 ) 
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Z12=imread( 'P12_0V.TIF' ) 
Z12=double (Z12) 

[Y12,X] =max (Z12) 
nunxl2=max  ( Y12 ) 
Z13=imread( 'P13_0V.TIF' ) 
Z13=double (Z13) 

[Y13,X] =max (Z13) 
nunxl3=max  (Y13) 
Z14=imread( 'P14_0V.TIF' ) 
Z14=double  (Z14) 

[Y14,X] =max (Z14) 
nunxl4=max  (Y14 ) 
Z15=imread( 'P15_0V.TIF' ) 
Z15=double (Z15) 

[Y15,  X] =max (Z15) 
nunxl5=max  (Y15) 
Z16=imread( 'P16_0V.TIF' ) 
Z16=double  (Z16) 

[Y16, X] =max (Z16) 
nunxl6=max  (Y16) 
Z17=imread( 'P17_0V.TIF' ) 
Z17=double  (Z17) 

[ Y17 , X] =max ( Z17 ) 
nunxl7=max  (Y17 ) 
Z18=imread( 'P18_0V.TIF' ) 
Z18=double  (Z18) 

[Y18,X] =max (Z18) 
nunxl8=max  (Y18 ) 
Z19=imread( 'P19_0V.TIF' ) 
Z19=double  (Z19) 

[Y19,X] =max (Z19) 
nunxl9=max  (Y19) 
Z20=imread( 'P20_0V.TIF' ) 
Z20=double  (Z20) 

[ Y20 , X] =max ( Z20 ) 
nunx2  0=max  (Y2  0 ) 
Z21=imread( 'P21_0V.TIF' ) 
Z21=double (Z21) 

[Y21,X] =max (Z21) 
nunx2  l=max  ( Y2 1 ) 
Z22=imread( 'P22_0V.TIF' ) 
Z22=double (Z22) 

[  Y22 , X] =max ( Z22 ) 
nunx22=max  ( Y22 ) 
Z23=imread( 'P23_0V.TIF' ) 
Z23=double (Z23) 

[Y23, X] =max (Z23) 
nunx23=max  (Y23) 
Z24=imread( 'P24_0V.TIF' ) 
Z24=double  (Z24) 

[Y24,X] =max (Z24) 
nunx2  4=max  (Y2  4 ) 
Z25=imread( 'P25_0V.TIF' ) 
Z25=double (Z25) 

[Y25, X] =max (Z25) 
nunx25=max  (Y25) 
Z26=imread( 'P26_0V.TIF' ) 


116 


Z26=double  (Z26) 

[Y26, X] =max (Z26) 
nunx2  6=max  (Y26) 

Z27=imread( 'P27_0V.TIF' ) 

Z27=double  (Z27) 

[ Y27 , X] =max ( Z27 ) 
nunx2  7=max  (Y2  7 ) 

Z28=imread( 'P28_0V.TIF' ) 

Z28=double  (Z28) 

[ Y28 , X] =max ( Z28 ) 
nunx2  8=max  (Y2  8 ) 

Z28=imread( 'P28_0V.TIF' ) 

Z28=double  (Z28) 

[ Y28 , X] =max ( Z28 ) 
nunx2  8=max  (Y2  8 ) 

Z29=imread( 'P29_0V.TIF' ) 

Z29=double  (Z29) 

[Y29, X] =max (Z29) 
nunx2  9=max  (Y29) 

Z30=imread( 'P30_0V.TIF' ) 

Z30=double  (Z30) 

[ Y30 , X] =max ( Z30 ) 
nunx30=max  ( Y30 ) 

Z31=imread( 'P31_0V.TIF' ) 

Z31=double  (Z31) 

[Y31,X] =max (Z31) 
nunx31=max  ( Y31 ) 

%POSITION  1-FROM  OV  TO  15V  (0.5V  INCREMENTS)  (DC) 
%READ  TIF  FILES  AND  EXTRACT  THE  MAX  VALUE 

Al=imread ( 'PI . TIF ' ) 

Al=double (Al) 

[B1 , X] =max (Al ) 
nunyl=max  (B1 ) 

A2=imread ( 'P2 . TIF  '  ) 

A2=double (A2) 

[B2 ,  X] =max  (A2 ) 
nuny2=max  (B2 ) 

A3=imread ( 'P3 . TIF ' ) 

A3=double (A3 ) 

[B3, X] =max (A3) 
nuny3=max  (B3 ) 

A4=imread ( 'P4 . TIF ' ) 

A4=double (A4) 

[B4 , X] =max (A4 ) 
nuny4=max  (B4 ) 

A5=imread ( ' P5 . TIF  '  ) 

A5=double (A5 ) 

[B5, X] =max (A5) 
mmy5=max (B5 ) 

A6=imread ( ' P6 . TIF ' ) 

A6=double (A6 ) 

[B6, X] =max (A6) 
nuny6=max  (B6 ) 

A7=imread ( ' P7 . TIF ' ) 
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A7=double (A7) 

[B7 , X] =max (A7 ) 
nuny7=max  (B7 ) 

A8=imread ( 'P8 . TIF ' ) 
A8=double (A8 ) 

[B8 , X] =max (A8 ) 
nuny8=max  (B8 ) 

A9=imread ( 'P9 . TIF ' ) 
A9=double (A9 ) 

[B9, X] =max (A9) 
nuny9=max  (B9 ) 
A10=imread( 'PIO.TIF'  ) 
A10=double (AlO) 
[B10,X] =max (AlO) 
nunylO=max  (BIO ) 
All=imread( 'Pll.TIF' ) 
All=double (All) 
[B11,X] =max (All) 
nunyl l=max  (B1 1 ) 
A12=imread( 'P12.TIF' ) 
A12=double (A12) 
[B12,X] =max (A12) 
mmyl2=max (B12 ) 
A13=imread( 'P13.TIF' ) 
A13=double (A13) 
[B13,X] =max (A13) 
nunyl3=max  (B13) 
A14=imread( 'P14.TIF'  ) 
A14=double (A14) 
[B14,X] =max (A14) 
nunyl4=max  (B14 ) 
A15=imread ( ' P15 . TIF ' ) 
A15=double (A15) 

[B15, X] =max (A15) 
nunyl5=max  (B15) 
A16=imread( 'P16.TIF'  ) 
A16=double (A16) 

[B16, X] =max (A16) 
nunyl6=max  (B16) 
A17=imread( 'P17.TIF'  ) 
A17=double (A17) 

[B17 ,  X] =max (A17 ) 
nunyl7=max  (B17 ) 
A18=imread( 'P18.TIF'  ) 
A18=double (A18) 

[B18 ,  X] =max (A18 ) 
nunyl8=max  (B18 ) 
A19=imread( 'P19.TIF'  ) 
A19=double (A19) 
[B19,X] =max (A19) 
nunyl9=max  (B19) 

A2  0=imread( 'P20.TIF'  ) 
A20=double (A20) 

[B2  0 ,  X] =max (A2  0 ) 
nuny2  0=max  (B2  0 ) 
A21=imread( 'P21.TIF'  ) 
A21=double (A21) 
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[B21,X] =max (A21) 
nuny2 l=max  (B2 1 ) 
A22=imread ( ' P22 . TIF '  ) 
A22=double (A22) 

[B22 , X] =max (A22 ) 
nuny22=max  (B22 ) 
A23=imread( 'P23.TIF'  ) 
A23=double (A23) 

[B23, X] =max (A23) 
nuny23=max  (B23) 
A24=imread( 'P24.TIF'  ) 
A24=double (A24) 

[B2  4 ,  X] =max (A2  4 ) 
nuny2  4=max  (B2  4 ) 
A25=imread ( ' P25 . TIF ' ) 
A25=double (A25) 

[B25, X] =max (A25) 
nuny25=max  (B25) 

A2  6=imread( 'P26.TIF'  ) 
A2  6=double (A2  6 ) 

[B26, X] =max (A26) 
nuny2  6=max  (B2  6) 
A27=imread ( ' P27 . TIF ' ) 
A27=double (A27) 

[B27 , X] =max (A27 ) 
nuny2  7=max  (B2  7 ) 

A2  8=imread( 'P28.TIF'  ) 
A28=double (A28) 

[B2  8 ,  X] =max (A2  8 ) 
nuny2  8=max  (B2  8 ) 

A2  9=imread( 'P29.TIF'  ) 
A2  9=double (A2  9 ) 

[B29, X] =max (A29) 
nuny2  9=max  (B2  9) 
A30=imread( 'P30.TIF'  ) 
A30=double (A30) 

[B30 ,  X] =max (A30 ) 
nuny30=max  (B30 ) 
A3I=imread( 'P31.TIF'  ) 
A3I=double (A3I) 
[B3I,X] =max (A3I) 
mmy3I=max (B3I ) 


%CREATE  MATRIX  FOR  MAX  VALUES  OF  0  VOLTS 

X=[mmxl  mmx2  mmx3  mmx4  mmx5  mmx6  mmx7  mmx8  mmx9  mmxIO  mmxll  mmxI2 
nunxI4  mmxIS  mmxI6  mmxI7  mmxI8  mmxI9  mmx20  mmx2I  mmx22  mmx23  mmx24 
nunx26  mmx27  mmx28  mmx29  mmx30  mmx3I] 


%CREATE  MATRIX  FOR  MAX  VALUES  OF  15  VOLTS 

Y=[mmyl  mmy2  mmy3  mmy4  mmy5  mmy6  mmy7  mmy8  mmy9  mmylO  mmyll  mmyI2 
nunyI4  minyI5  mmyI6  mmyI7  mmyI8  mmyI9  mmy20  mmy2I  mmy22  mmy23  mmy24 
nuny2  6  mmy2  7  mmy2  8  mmy2  9  mmy30  mmy3I] 


mmxI3 
mmx2  5 


mmyI3 
mmy2  5 


%CALCULATE  NORMALIZED  INTENSITY 
Io=240 

1= (Y-Io) . / (X-Io) 
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%SET  VALUES  FOR  APPLIED  VOLTAGE 
V=0:0.5:15 


%EVALUATE  APPLIED  VOLTAGE  W.R.T.  DISTANCE  BETWEEN  PLATES  (V/pm) 
%SET  VALUES  FOR  OUTPUT  CURRENT 


"6 

Iout= 

=  [0. 

001 

0 

.032  0 

.  537 

1 . 

049 

4 . 080 

6.335  8 

.  710 

11 

.  158 

13 

.  658 

16, 

.200 

18 

.767 

21 , 

.355 

2 

3 

.  952 

26.5 

62 

29. 

190 

31 

.836  34. 

503 

37  . 

187 

39. 

.895 

42  , 

.  647 

45 

.  420 

48  , 

.390 

5 

1 

.230 

54 . 0 

90 

56. 

960 

59 

.890  62. 

760 

65. 

640 

68  . 

.540 

71 , 

.340 

74 .230] 

%  P=V.*Iout 
%  R= (V. *V) . /P 
%  VI=Iout . *R/I9I 
VI=V/I9I 


%CREATE  THE  GRAPH 
plot (V, I , ' bo ' ) 
grid  on 

xlabel (' Applied  Voltage  (V) ') 
ylabel (' Intensity  Ratio') 

title (' Normalized  Intensity  Ratio  vs  Applied  Voltage  (Experimental 
Data) ' ) 


%EVALUATION  OF  THEORITICAL 
CORRESPONDING  GRAPH 


MODEL  FOR  DC  AND  CREATION  OF  THE 


[0 

0 . 001 

0 . 002 

0 . 003 

0 . 004 

0 . 0044 

0 . 005 

0.006 

0 . 007 

008 

0 .01 

0 .012 

0 . 0132 

0 . 015 

0.017 

0 . 02 

0 . 021 

0 . 023 

025 

0 . 027 

0.029 

0 . 0314 

0 . 033 

0 . 035 

0 . 037 

0 . 039 

0.0412 

043 

0 . 045 

0 . 047 

0 . 049 

0 . 0521 

0 . 054 

0 . 056 

0 . 058 

0.06]; 

Iold= [0 . 0173931861141514 
0 .0171766638980226 
0 . 0168149071852186 
0 .0162940163629757 
0 . 0150242440677031 
0 .0131890887893245 
0 .0121406835833828 
0 . 0109547598726452 
0 .0100454544116619 
0 . 00921432745937045 
0 . 00850034248019468 


0 . 0173686237268924 
0 .0170149206611566 
0 .0166530028607825 
0 .0158421498452957 
0 .0145185911735649 
0 .0129837833235535 
0 .0117481446599407 
0 .0106754611667708 
0 . 00975404133410867 
0 . 00896439018545242 


0 . 00818920174584799 
0.00783888525452678  0.00767012138997172  0.00750822234977514 
Inorm=Iold/0 .01739318611415140 


0 . 0172957143671621 
0 .0169392356157452 
0 . 0164875036781727 
0 .0153323559252167 
0 . 0139712700358284 
0 .0125496590763713 
0 .0113730073666232 
0 . 0103523335916252 
0 . 00945037407473868 
0 . 00872663147449506 
0 . 00801488318079925 


n=2 ; 

p=polyf it (Inorm, E, n) 
xi= (0 . 40:0.01:1.0); 
yi=polyval (p, xi) ; 


figure 

plot ( Inorm, E, 'o',xi,yi, '-') 
grid  on 

xlabel (' Intensity  Ratio '), ylabel (' E-Field  (V/mum') 
title ('Curve  Fit  for  DC  Theoretical  Model') 
legend ( ' Theoretical ' ,  ' Curve  Fit '  ) 

text  (0.42,0.005,  'Etheory=  (0 . 0842*I''2)-(0. 21 62*1) +0.1351') 
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%APPLY  CORRECTION 

Etheory= (0 . 0842*1. *I)-(0. 21 62*1) +0.1351 

%CREATE  THE  FINAL  GRAPH 
figure 

plot (V, Etheory, 'o') 
hold  on 
plot (V, VI, ' - ' ) 
hold  off 

legend ( ' Experimental ' , ' Theoretical ' ) 
grid  on 

xlabel (' Applied  Voltage  (V)  '),  ylabel (' E-Field  (V/\mum)') 
title (' E-Field  vs  Applied  Voltage') 

%CREATE  THE  FINAL  GRAPH  (Log  Scale) 
figure 

semilogy (V, Etheory,  '-o') 
hold  on 

semilogy (V, VI , '-') 
hold  off 

legend ( ' Experimental ' , ' Theoretical ' ) 
grid  on 

xlabel (' Applied  Voltage  (V)  '),  ylabel (' E-Field  (V/\mum)') 
title (' E-Field  vs  Applied  Voltage') 

%CREATE  THE  FINAL  GRAPH  (Log  Scale) 
figure 

loglog (V, Etheory, '-o') 
hold  on 

loglog (V, VI , ' - ' ) 
hold  off 

legend ( ' Experimental ' , ' Theoretical ' ) 
grid  on 

xlabel (' Applied  Voltage  (V) '), ylabel (' E-Field  (V/\mum)') 
title (' E-Field  vs  Applied  Voltage') 

%CREATE  THE  FINAL  GRAPH  WITHOUT  LINE 
figure 

plot (V, Etheory, '-o') 
legend ( ' Experimental ' ) 
grid  on 

xlabel (' Applied  Voltage  (V) '), ylabel (' E-Field  (V/\mum)') 
title (' E-Field  vs  Applied  Voltage') 
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APPENDIX  D.  MATLAB  CODE  FOR  COMSOL  DATA 


%ANALYZE  DATA  FROM  COMSOL 
%CREATED  BY  PAVLOS  ANDRIKOPOULOS 

%READ  DATA  FROM  .txt  FILES 

load  EX. txt 

load  EY.txt 

load  EZ.txt 

xl=EX  (  : , 2) 

yl=EX(:,l) 

x2=EY ( : , 2) 

y2=EY(:,l) 

x3=EZ  (  : , 2) 

y3=EZ ( : , 1) 

x=sqrt (xl . *xl+x2 . *x2+x3 . *x3) 
d=yl 

%CREATE  THE  GRAPHS 

%CREATE  GRAPH  FOR  ELECTRIC  FIELD  IN  X-DIRECTION 

figure 

plot (yl , xl ) 

grid  on 

xlabel (' Distance  (m) ') 

ylabel (' Electric  Field  (V/m) ') 

title  (' Electric  Field  -  x  component') 

%CREATE  GRAPH  FOR  ELECTRIC  FIELD  IN  Y-DIRECTION 

figure 

plot (y2 , x2 ) 

grid  on 

xlabel (' Distance  (m) ') 

ylabel (' Electric  Field  (V/m) ') 

title  (' Electric  Field  -  y  component') 

%CREATE  GRAPH  FOR  ELECTRIC  FIELD  IN  Z-DIRECTION 

figure 

plot (y3, x3) 

grid  on 

xlabel (' Distance  (m) ') 

ylabel (' Electric  Field  (V/m) ') 

title  (' Electric  Field  -  z  component') 

%CREATE  GRAPH  FOR  MAGNITUDE  OF  TOTAL  ELECTRIC  FIELD 

figure 

plot (d, x) 

grid  on 

xlabel (' Distance  (m) ') 
ylabel (' Electric  Field  (V/m) ') 
title (' Electric  Field  (V/m) ') 
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